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^ ! We investigate how the empirical properties of hot X-ray-emitting gas in a 

^ ! sample of seven starburst and three normal edge-on spiral galaxies (a sample 

Q\ ! which covers the full range of star-formation intensity found in disk galaxies) 

^ ! correlate with the size, mass, star formation rate and star formation intensity in 

O ! the host galaxies. From this analysis we investigate various aspects of mechanical 

2 [ energy "feedback" — the return of energy to the ISM from massive star super- 

r-| ! novae and stellar winds — on galactic scales. The X-ray observations make use of 

the unprecedented spatial resolution of the Chandra X-ray observatory to remove 
O I X-ray emission from point sources more accurately than in any previous study, 

^ [ and hence obtain the X-ray properties of the diffuse thermal emission alone. In- 

I triguingly, the diffuse X-ray properties of the normal spirals (both in their disks 

. ^ I and halos) fall where extrapolation of the trends from the starburst galaxies with 

^ I superwinds would predict. We demonstrate, using a variety of multi- wavelength 

\ star formation rate and intensity indicators, that the luminosity of diffuse X-ray 

emission in the disk (and where detected, in the halo) is directly proportional to 
the rate of mechanical energy feedback from massive stars in the host galaxies. 
Accretion of gas from the IGM does not appear to be a significant contributor 
to the diffuse X-ray emission in this sample. Nevertheless, with only three non- 
starburst normal spiral galaxies it is hard to exclude an accretion-based origin for 
extra-planar diffuse X-ray emission around normal star-forming galaxies. Larger 
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galaxies tend to have more extended X-ray-emitting halos, but galaxy mass ap- 
pears to play no role in determining the properties of the disk or extra-planar 
X-ray emitting plasma. The combination of these luminosity and size correla- 
tions leads to a correlation between the surface brightness of the diffuse X-ray 
emission and the mean star formation rate per unit area in the disk (calculated 
from the Far- Infrared luminosity and the optical size of the galaxy, Lfir/-D25)- 
Further observational work of this form will allow empirical constraints to be 
made on the critical star formation rate per unit disk area necessary to blow hot 
gas out of the disk into the halo (Mac Low & McCray 1988). We show that a 
minor generalization of standard superbubble theory directly predicts a critical 
star formation rate per unit area for superbubble blow out from the disk, and by 
extension for superwinds to blow out of the gaseous halos of their host galaxy. 
At present there are a variety of poorly-known parameters in this theory which 
complicate comparison between observation and theory, making it impossible to 
assess the quantitative accuracy of standard superbubble blow out theory. We 
argue that the crucial spatial region around a galaxy that controls whether gas in 
starburst-driven superwinds will escape into the IGM is not the outer halo ~ 100 
kpc from the host galaxy, but the inner few halo scale heights, within ~ 20 kpc 
of the galaxy plane. Given the properties of the gaseous halos we observe, super- 
wind outflows from disk galaxies of mass M ~ 10^° - 10^^ M© should still eject 
some fraction of their material into the IGM. 

Subject headings: ISM: jets and outflows — ISM: bubbles — galaxies: individual 
(NGC 253; NGC 891; NGC 1482; NGC 3034 (M82); NGC 3073; NGC 3079; 
NGC 3628; NGC 4244; NGC 4631; NGC 4945; NGC 6503) — galaxies: halos — 
galaxies: starburst — X-rays: galaxies 

1. Introduction 

Massive stars exercise a profound influence over the baryonic component of the Universe, 
through the return of ionizing radiation, and via supernovae (SNe), kinetic energy and metal- 
enriched gas, back into the interstellar medium (ISM) from which they form — usually called 
"feedback" . Feedback influences gas- phase conditions in the immediate environment of the 
clusters within which the massive stars form (McKee 1995; HoUenbach & Tielens 1997; 
Wiseman & Ho 1998; Pudritz & Fiege 2000), on galactic scales the phase structure and 
energetics of the ISM (McKee & Ostriker 1977; Cox 1981; Norman & Ikeuchi 1989; Norman 
& Ferrara 1996), and on multi-Mpc scales the thermodynamics and enrichment of the inter- 
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galactic medium (IGM — e.g. see Chiang, Ryu & Vishniac 1988; Heckman, Armus Sz Miley 
1990; Shapiro, Giroux & Babul 1994; Voit 1996; Heckman 1999; Aguirre et al. 2001). 

The vast range of spatial scales involved is only one of the difficulties encountered 
in attempting to study feedback. Restricting the discussion to purely mechanical feedback 
from SNe and stellar winds (commonly termed SN feedback) , a further difficulty is the broad 
range of complicated gas-phase physics - (magneto)hydrodynamic effects such as shocks and 
turbulence, thermal conduction, and non-ionization equilibrium emission processes. A final 
complication is that much of the energy and metal-enriched material involved is in the hard- 
to-observe coronal (T > 10^ K) and hot (T > 10^ K) gas phases, necessitating the use of 
space-based FUV and X-ray telescopes. 

X-ray observations of the halo-regions of nearby edge-on disk galaxies provide the best 
single probe of the action of mechanical feedback on the galactic scale. Optical spectroscopy 
and imaging, sensitive to warm neutral and ionized gas with T ~ 10'^ K, are a less direct probe 
of mechanical feedback as much of the excitation energy is supplied by ionizing radiation 
from massive stars (radiative feedback), even when the kinematics and spatial distribution 
of this gas are controlled by SN feedback. If mechanical feedback processes transport mass, 
newly-synthesized metals, and energy from galaxies into the IGM, then this material must 
pass through the halo of the host galaxy on its way into the IGM, where its properties can 
be ascertained before it expands and fades into observational invisibility. 

Diffuse thermal X-ray emission is seen to extend to large distances (5 to 20 kpc) above 
the planes of edge-on starburst galaxies (e.g. Fabbiano, Heckman & Keel 1990; Armus 
et al. 1995; Read, Ponman & Strickland 1997; Dahlem, Weaver & Heckman 1998), which 
also show unambiguous optical emission line evidence for energetic galactic scale outfiows 
— "superwinds" (Heckman et al. 1990). These outflows are potentially very signiflcant in 
enriching and heating the IGM (Heckman 1999). 

Both fully-numerical and semi-analytical theoretical modeling (Suchkov et al. 1994; 
Silich & Tenorio-Tagle 1998; Murakami & Babul 1999) suggest that the properties of the 
gaseous halos these superwinds expand into play a major role in determining the fate of the 
material in superwinds (i.e. escape into the IGM, or confinement and eventual return to the 
disk). Although simulations of starburst-driven winds should be treated with caution, given 
the lack of calibration against observational data, one robust result is the sensitivity of the 
wind properties (in particular mass and metal loss to the IGM) on the disk and halo ISM 
properties (e.g. De Young & Heckman 1994; Suchkov et al. 1994; Sihch & Tenorio-Tagle 
1998; Strickland & Stevens 2000; Sihch & Tenorio-Tagle 2001). In our view, the emphasis 
given to galactic mass as a primary parameter infiuencing gas and metal ejection efficiencies 
(e.g. Dekel & Silk 1986; Mac Low & Ferrara 1999; Ferrara & Tolstoy 2000; Ferrara et al. 
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2000) blinds non-expert readers to the more fundamental role played by the poorly-known 
disk and halo ISM distributions. The gravitational potential only enters at a secondary level, 
by influencing the gas distribution and flnally via the escape velocity. Redressing the lack of 
observational knowledge about the gaseous halos that superwinds must pass through is one 
motivation for studying the halos of normal star-forming galaxies. 

Modern theories of the ISM in normal spiral galaxies also predict hot gas in the halo, 
due to the interchange of material between the disk and halo through "galactic fountains" 
(Shapiro & Field 1976; Bregman 1980; de Avillez 2000), or in a more-localized form through 
"chimneys" (Norman & Ikeuchi 1989). Energy feed back into the ISM from supcrnovac is 
believed to create bubbles of hot, metal-enriched, gas, which blow out of the disk and vent 
their contents into the halo of the galaxy. After > 10^ years or so, this material cools and 
falls back to the disk, possibly in a form analogous to the Galactic high velocity clouds. 

Unfortunately, many of the fundamental aspects of these fountain/ chimney disk/halo 
interaction theories are observationally unconstrained, fundamentally due to our current lack 
of knowledge about the X-ray properties of spiral galaxy halos. In the theoretical paradigm, 
fountains are primarily driven by the over-pressurized hot gas created by SNe, and the most 
direct diagnostic of the energetics, and elemental composition, of this gas is the thermal 
X-ray emission it produces. 

Studies of nearby edge-on star-forming galaxies with Einstein and the ROSAT PSPC 
revealed only one normal spiral, NGC 891 (Bregman & Pildis 1994), with X-ray-emitting 
gas in its halo^. Although edge-on normal spiral galaxies did not receive as much attention 
with Einstein & ROSAT as other classes of galaxy, those observations that were performed 
failed to produce believable detections of hot halo gas in other edge-on^ normal star-forming 
galaxies (Bregman & Glassgold 1982; Vogler, Pietsch & Kahabka 1995; Bregman & Houck 
1997; Fabbiano & Juda 1997; Benson ct al. 2000). X-ray studies of our own galaxy's halo 
and disk/halo interface (e.g. Kuntz & Snowden 2000) provide information complementary 



^NGC 4631 is often presented as an example of a "normal" spiral galaxy with a X-ray-emitting halo (e.g. 
see Wang et al. 2001). However, its large-scale- height non- thermal radio halo, warm IRAS 60 to 100/um flux 
ratio, high Ha luminosity and optically disturbed morphology show that it is a highly atypical spiral galaxy. 
It is a good example of a galaxy experiencing a mild, disk- wide, starburst event (GoUa & Wielebinski 1994). 
We believe that the ROSAT PSPC-based detection of a = 2.6 x 10^9 diffuse X-ray halo around NGC 
4565, claimed by Vogler et al. (1995), is incorrect, as the archival 2.9 ksec Chandra ACIS-S observation 
(ObsID 404) reveals that the X-ray emission ^ 2' cast of the nucleus (that comprised the majority of the 
emission originally considered to be diffuse) is due to two point sources. 

^Diffuse X-ray emission in approximately face-on star-forming galaxies (e.g. Snowden & Pietsch 1995; 
Cui et al. 1996; Ehle et al. 1998), can not be unambiguously associated with a hot halo, for obvious reasons. 
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to the study of other galaxies, but are difficult to perform and interpret due to our location 
within the disk. 

More progress on investigating the disk/halo interface in normal galaxies has been made 
at other wavelengths, primarily in the optical, where extra-planar Ha emission (e.g. Dettmar 
1992; Rand 1996; Hoopes, Walterbos & Rand 1999; Collins et al. 2000; Rossa & Dettmar 
2003; Miller & Veilleux 2003) and dust (Howk & Savage 1999) has been detected in a 

variety of nearby, normal spiral galaxies. These studies have shown that the presence of such 
tracers of disk/halo interaction increases with higher levels of star-formation intensity in the 
underlying disk, but the mechanisms behind these phenomena, and their energetics, are still 
unclear. 

It is also not clear to what extent fountains are a distinct phenomenon from superwinds. 
It is not unreasonable to think that physical conditions and kinematics of the gas making up 
a fountain may be quite different from that in a superwinds, despite the general similarity in 
driving mechanism. For example, the filling factor of the dominant X-ray-emitting plasma 
might be high in a fountain (Bregman 1980), as compared to superwinds where the filling 
factor is low (Chevalier & Clegg 1985; Suchkov et al. 1994; Strickland & Stevens 2000; 
Strickland et al. 2000). Another intriguingly hint of a difference comes from the ROSAT 
PSPC-based sample of both normal spiral galaxies and starburst galaxies of all inclinations 
by Read & Ponman (2001). They find that the X-ray-emitting gas "in the coronal systems 
[i.e. galaxies assumed to have a hydrostatic hot gas corona/halo] appears to be cooler than 
that seen in the wind systems [i.e. starburst galaxies with superwinds]". 

Practitioners of semi- analytical (e.g. Benson et al. 2000), and numerical cosmological 
models (e.g. Toft et al. 2002) of galaxy formation and evolution offer a radically different 
origin for X-ray emitting halos around spiral galaxies. Inter-galactic gas continually flows 
into galactic halos, is heated in an accretion shock to the halo virial temperature (although 
Katz et al. (2003) discuss simulations in which only a small fraction of the accreted gas 
is heated to halo virial temperatures), and then cools and accretes onto the disk. These 
models predict red-shift zero spiral galaxy X-ray luminosities at, or just below, the ROSAT- 
based observational upper limits. The predicted X-ray luminosity is a strong function of the 
galaxy mass, Lx oc v^^^, and that the temperature of the gas is similar to the halo virial 
temperature kT ~ 0.5 /imn I'rot- The majority of the diffuse X-ray emission arises close to 
the disk, the X-ray surface brightness dropping off semi-exponentially with height above the 
disk. These halo properties should be largely independent of the host disk's SF intensity, 
and the exact numerical implementation of SN feedback (J. Sommer-Larsen, 2002, private 
communication) . 



Table 1. Basic physical properties of the sample galaxies 
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Note. — Column 2: Coordinates of the dynamical center of the galaxy, where available. Column 3: Inclination. Column 4: Position angle of the galactic disk. Multiple values 
for the inclination and position angle are given when several values are often reported in the literature. Column 5: Inclination-corrected circular velocity within the disk. ^1 umn 
6: Assumed distance. Unless a specific reference is given, distances were calculated using the galaxy recessional velocity with respect to the microwave background (de Vaucouleurs 

et al. 1991), and Hq = 75kms^-'^ Mpe^^. Column 7: Physical distance corresponding to an angular size of 1 areseeond. Column 8: IRAS 60/jrri flux in Janskys. Column 9: IRAS 
60 to 100/im flux ratio. Column 10: Total IR luminosity based on the observed IRAS fluxes, Lju = 5.67 X lO^-D^p^. X (13.48/i2 + 5.I6/25 + 2.58/60 + /lOo) Lq) , where Dmpc is the 
distance to the galaxy in Mpc (Sanders & Mirabel 1996). Given the large angular size of these local galaxies we use IRAS fluxes from extended source analyses, the exact source 
of the IRAS data used is given in a case by case basis above. Where multiple measurements were available for any galaxy, we used the following sources in order of preference: 
Rice et al. (1988), and Soifer et al. (1989). Column 11: K-band luminosity, calculated from the K-band total magnitudes given in Jarret et al. (2003) and converted to the 
Cousins-Glass- Johnson system using the transformations presented in Carpenter (2001). Zero points for conversion from magnitude to flux were taken from Bessel, Castelli & Plez 
(1998). Note that the units are Lq, not Lk 0. Column 12: B-band luminosity from RC3 catalog Bx magnitude (de Vaucouleurs et al. 1991), and corrected for Galactic extinction 
using reddening values from Burstein & Heiles (1982). Note that no correction for extinction internal to the speciflc galaxy has been made. Column 13: Observed Ha-luminosity 
in units of 10*" ergs~^ . Note these values have not been corrected for extinction. Measurements that have not been corrected for [Nil] emission are marked with a preceding A'^. 
Column 14: Baryonic mass (stellar plus Hi) estimated from the K-band TuUy Fisher relationship of Bell & de Jong (2001), where Mtf = lO^ '''^ x (i^rot/lOOkms"! )3 Si Mg. 

References. — (a) Weliachew, Fomalont & Greisen (1984), (b) Ichikawa et al. (1995), (c) Lynds & Sandage (1963), (d) Lehnert & Heckman (1995), (e) Preedman et al. 
(1994), (f) Soifer et al. (1989), (g) McCarthy, Heckman & van Bruegel (1987), (h) 2MASS Second Incremental Release, (i) Hameed & Devereux (1999), (k) Calculated assuming 
v^ot = W2o/2sin2 j, using the W20 value from Roth, Mould & Davies (1991), (1) Hameed & Devereux (1999), (m) Sorai et al. (2000), (n) Pence (1981), (o) Puche, Carignan 
& van Gorkom (1991), (p) Puche & Carignan (1988), (q) Prom Ha images published in Strickland et al. (2002), (r) Douglas et al. (1996), (s) Irwin & Sofue (1996), (t) Sofue 
(1997), (u) Soifer et al. (1987) (v) Fabbiano et al. (1990), (w) Irwin & Seaquist (1991), (x) Ott et al. (2001), (y) Koorneef (1993), (z) Rice et al. (1988), (aa) Dynamical center, 
GoUa & Wielebinski (1994), (ab) Rand (1994), (ac) GoUa (1999), (ad) Hoopes et al. (1999), (ae) Central optical & X-ray source. Lira, Johnson & Lawrence (2002), (af) de 
Vaucouleurs et al. (1991), (ag) Bottema & Gerritsen (1997), (all) Karachentsev &: Sharina (1997), (ai) Kennicutt (1998a). (aj) Nuclear 1.4 GHz continuum point source, Rupen 



(1991), (ak) Hi dynamical center, Oiling (1996). 
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In order to investigate some of these questions we have compiled a sample of 10 edge-on 
star-forming galaxies which have both ground-based narrow-band optical imaging (probing 
extra-planar gas at T ~ 10^) and X-ray observations with the Chandra X-ray Observatory. 
The spatial resolution of the Chandra telescope (~ (X'S FWHM) is an order of magnitude 
better than any previous or currently-existing X-ray instrument. As approximately half 
of the X-ray luminosity of a star-forming galaxies is due to X-ray luminous point sources, 
accurate studies of the X-ray diffuse emission require high spatial resolution in order to allow 
point source removal. 

In Strickland ct al. (2004, henceforth referred to as Paper 1) we described the sample, 
the X-ray and optical data reduction performed, and presented results and interpretation of 
a spatial and spectral analysis of the diffuse X-ray emission. 

In this paper we investigate how the empirical properties of the hot gas, in particular 
gas in the halos of these galaxies, correlate with the size, mass, star formation rate and star 
formation intensity in the host galaxies. Prom this analysis we investigate various aspects 
of mechanical energy "feedback" — the return of energy to the ISM from massive star 
supernovae and stellar winds — on galactic scales. 

We describe the sample and their star-formation properties in § 2. In § 3 we provide 
a brief summary of the results from Paper I, before moving onto considering what correla- 
tions exist between the properties of the diffuse X-ray emission and other multi- wavelength 
properties of the galaxies in § 4. In § 5 we discuss these results with respect to various 
aspects of feedback-related theory, in particular the conditions necessary for blow-out of gas 
from the disk and halo in normal galaxies and in starbursts, and the fraction of supernova 
mechanical energy available to drive large-scale ISM motions. Our results and conclusions 
are summarized in § 6. 

2. Sample galeixies 

The sample is formed of seven starburst galaxies and three "normal" star-forming galax- 
ies, all approximately edge-on {i > 60°) moderate luminosity (Lbol ~ 10^° L©)) moderate- 
mass (M ~ 10^*^ - 10^^ Mq) disk galaxies within a distance of D < 20Mpc. This sample 
contains all the disk galaxies for which Chandra ACIS data was available to us as of 2002 
September, including both public data and our own propriety observations. Although the 
sample size is small, in particular with respect to non-starburst systems, it does span the 
full range of star formation activity found in disk galaxies. The starburst galaxies in our 
sample are all good examples of the "typical" starburst galaxy in the local universe, i.e. 
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infrared warm galaxies with a far-infrared luminosity LpiR ~ a few x 10^° Lq (Soifer et al. 
1987). The three non-starburst systems in the sample are NGC 891, NGC 6503 and NGC 
4244. NGC 891 is a well-studied spiral galaxy, believed to be similar to our own galaxy in 
mass and star formation rate. It is the only "normal" spiral galaxy (apart from possibly our 
own) for which robust evidence of extra-planar diffuse X-ray emission has been discovered 
(Bregman & Pildis 1994). Both NGC 6503 and NGC 4244 are relatively low mass spiral 
galaxies. We refer the reader to Paper I for a more detailed discussion of the sample and the 
distinction between starburst and normal spiral galaxies. 

Basic physical properties of the sample galaxies are shown in Table 1. We distinguish 
between starburst galaxies and more-normal star-forming galaxies on the basis of their IR- 
warmth — specifically the IRAS 60 to 100//m flux ratio /eo/Zioo- Higher star formation 
rates per unit disk area lead to higher dust temperatures, and hence higher /eo/Zioo ratios, 
with starburst galaxies typically having /eo/Zioo > 0.4. We will refer to those galaxies with 
/eo/Zioo < 0.4 as the normal spirals. 

Intuitively, we would expect the effect of mechanical energy feedback on the galactic 
scale to depend on the total star formation rate, and/or the star formation rate per unit 
volume or area. Star formation and supernova rates, and rates per unit area, for the sample 
galaxies arc given in Table 2. Note that galaxies arc listed in all tables in order of decreasing 
/eo/Zioo ratio, i.e. approximately in order of decreasing star formation rate per unit area. 

2.1. Star formation rates and intensities 

We also consider two other simple, and widely-available proxies of the star-formation 
intensity: FIR and non-thermal radio face-on surface brightnesses, quantifled for our sam- 
ple in Table 2. Although somewhat simplistic estimators of SF intensity, these do have 
the advantages of being available for almost all local disk galaxies from uniformly-obtained 
datasets, and are distance-independent. 

Very detailed studies of the star-formation rates, spatial distributions, and history, are 
available for the brightest starbursts such as M82 (Ricke et al. 1993; McLeod et al. 1993; 
Satyapal et al. 1997; de Grijs, O'Connell & Gallagher 2001; Forster Schreiber et al. 2001) and 
NGC 253 (Engelbracht et al. 1998), but measurements of similar quality are not available 
for many of the starburst galaxies in our sample, let alone the normal spiral galaxies. We 
did calculate star formation intensities for all our sample galaxies based on all available 
optical, NIR and radio measurements we could flnd in the literature, but found the scatter 
in estimated SF intensity for any particular galaxy to be intolerably large. 



Table 2. Star- formation and supernova rates, SF intensity, and proxies thereof 
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0.0543 


1.9 X 


10-'' 


2.48 


7.70 


NGC 6503 


> 0.06 


0.20 


5.6 


8.4 


0.77 


1.18 


55.1 


0.492 




16.7 


0.0012 


41.9 


6.74 X 10^ 


34.6 


0.87 




0.0088 


4.0 X 


10^ 


2.97 


8.23 


NGC 891 


> 0.26 


4.2 


9.3 


26.1 


1.47 


4.10 


355.6 


1.14 




39.5 


0.049 


285.7 


1.50 X 10* 


55.2 


2.57 




0.0234 


1.9 X 


10^ 


3.19 


8.19 


NGC 4244 


> 0.03 


0.04 


10.2 


10.7 


1.85 


1.94 


27.2 


0.071 




3.5 


0.0004 


20.3 


5.65 X 10^ 


158.5 


2.84 




0.0003 


1.2 X 


lOl 


4.53 


7.72 


Note. — 


Column 2: 


Lower limits on 


the 


total galactic 


star formation rate 


( A^Q yr ^ ), based on 


the HcK luminosity and using the 


formulae given in Kennicutt (1998b). 


Note that no 


correction for extinction has been 


applied. 


Column 


3: Estimated 


total galact 


ic 


star format 


ion rate 


( Mo, yr- 


^), based on 


the IR lum 


inosity, 


where SFR/j^ — 


: 4.5 X 1(1 


l-**t 


IR [orgs^l ], 


again using 


the Kennicutt (1998b) formulae. Columns 4 


and 5: 


Inclination- 


-corrected diameter of the stellar disk in arcminutes (4) and kp 


c (5), based on 


the D25 


values given 


in the 


RC3 and the distances given in Table 1 


. The method of Tully & Fouqu« (1985) 


was used to correct the observed D25 values for inclination. 


Columns 6 and 7: K-band half light radii in 



arcminutes (6) and kpc (7). Original values were taken from the 2MASS Large Galaxy Atlas (Jarret et al. 2003). Column 8: IRAS Far-IR flux = 2.58 X /go + /lOO) Jy- To convert 

to units of ergs""*^ cm~^ multiply by 1.26 X 10""'^"'^. Column 9: Far-IR luminosity per unit stellar disk area /FIR/-D25, in distance-independent units of mjy arcsec"^. To convert to the 
commonly used units of lO***^ ergs""*^ kpc""*^ multiply /pir/^ 25 6.484. Column 10: One estimate of the mean galactic supernovae rate per unit area, -F'sN,FIR,D25 — ^SN / ^2^-' units 

of SN Myr""*" kpc~^, is based on the star-formation rate calculated from i^ip . Column 11: Estimated total galactic core-collapse supernova rate (yr""*"), where T^-sjv — O-^-C'/jf/lO"'""'" .Lq 
(Heckman et al. 1990). Column 12: 1.4 GHz radio flux in mJy, based on the NVSS data. See § 2.2 for more details. Column 13: 1.4 GHz luminosity in units of Solar luminosities, where 
Li 4GHz — 4:71 U"^ u f 1 4GHz ■ The major-axis half-light radius, in units of arcseconds and kpc, that we derive the NVSS data is given in columns 14 and 15 respectively. Column 16: Effective 
1.4 GHz surface brightness, in mJy arcsec"^. C^olumn 17: Another estimate of the mean SN rate per unit area -F^n pjpj 4CH ~ ^S^'/'^^l 4GHz (^S^i^i units of SN Myr""'" kpc~^), 
based on the star- formation rate calculated from but using NVSS-based size 61,4. Note that NGC 4945 was not observed in the NVSS. For NGC 4945 we use the nuclear region 1.4 

GHz flux and size from Elmouttie et al. (1997), whose observations had similar spatial resolution to that of the NVSS. Column 18: Logarithm of the gas consumption timescale in Myr, 
where log Tgas = log ^gas — log ^SFRj and Sgas and Sgpp_ are the surface density of gas ( Mq pc~^ ) and star- format ion rate ( Mq yr ~ ^ kpc~^ ) respectively. Data are from Kennicutt 
(1998a), except for NGC 1482, NGC 3628, NGC 4244, NGC 4631 and NGC 4945. For these galaxies the star formation rate was estimated from the IRAS fluxes. Gas masses were obtained 
from Elfhag et al. (1996) for NGC 1482 (using the CO brightness to H2 column density scaling in Kennicutt (1998a), and a radius of 2 kpc, as ISO observations by Dale et al. (2000) 
demonstrate that 80% of the star-formation occurs within this radius), Irwin & Sofue (1996) for NGC 3628 (central, R < 340 pc, H2 gas mass), from Oiling (1996) for NGC 4244 (total 
HI gas mass), from Golla & Wielebinski (1994) for NGC 4631 {H2 gas mass for R < 2.5 kpc) and from Dahlem et al. (1993) for NGC 4945 (molecular gas mass within central 1.1 kpc). 
Column 19: Logarithm of the mean mass surface density fi — Mq^-p / D2q- 
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The FIR luminosity divided by the optical isophotal diameter (in kpc) at 25th magni- 
tude, LfirZ-DIj, is often used as a proxy for the mean star-formation rate per unit disk area 
in spiral galaxies (we shall occasionally refer to the star formation rate per unit disk area 
as the SFRl for convenience). The presence of extra-planar Ha and/or radio emission, and 
extra-planar dust, is known to qualitatively correlate with this parameter (Rand 1996, see 
also Dettmar 1998; Howk & Savage 1999; Rossa & Dettmar 2000; Dahlem et al. 2001). 
It under-estimates the true SF intensity in starburst galaxies, as the star-formation is often 
much more concentrated than the old stellar hght measured by D2^. This method also over- 
estimates the mean SF intensity in galaxies with star-formation rates somewhat less than 
that of the Milky Way, where an increasing fraction of LpiR is due to cool dust heated by 
diffuse star-light, and not by massive stars. 

Non-thermal radio emission, due to synchrotron emission from relativistic electrons 
accelerated in supernova remnants, provides a good tracer of SF activity on galactic scales 
(Condon 1992). The 1.4 GHz luminosity of a galaxy is almost directly proportional to 
the SF rate, which, combined with the measured angular size of the radio-emitting region, 
provides a good SF intensity measurement (Dahlem et al. 2001) . 

The radio fluxes, sizes and SF intensity values given in Table 2 are based on 1.4 GHz 
images from the NRAO/VLA Sky Survey (NVSS, see Condon et al. 1998), with the exception 
of southern galaxy NGC 4945, which lies outside the region of sky surveyed by the NVSS. 
A description of how these values were obtained is given in § 2.2 

2.2. Radio fluxes, sizes and SF intensity measurements 

Radio fluxes and half-hght radii were measured from the processed NVSS 1.4 GHz 
radio images (http://www.cv.nrao.edu/nvss/postage.shtml). Fluxes were summed within 
rectangular apertures encompassing the optical disk of each galaxy, each of total width equal 
to the optical D25 diameter and breadth equal to the minimum of 1'5 (twice the angular 
resolution of the NVSS) or the angular equivalent of 4 kpc, and centered on the nuclear 
positions given in Table 1. The major-axis half-light radii 9ia were measured from the radio 
flux within the disk regions defined above using Ar = 5"-wide pixels. Inspection of the 
resulting radio surface brightness slices demonstrated that the fluxes at the extreme ends of 
the slices were in all cases statistically consistent with zero, allowing us to define the 100%- 
enclosed light radii as the edges of the surface brightness profiles. Given the linear, rather 
than radial, nature of these major-axis surface brightness profiles we obtained two half-light 
radii on either side of the nucleus for each galaxy, 9+r and 6'_r, which were combined to give 
the final half light radius ^1.4 = {{61, + Ol,)/ 2}^/'^ shown in Table 2. 
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We did not use the fluxes and sizes given in the NVSS Source Catalog itself, which 
quotes fluxes and sizes from a a 2- dimensional Gaussian model fit to the radio images. The 
NVSS Source Catalog underestimates the total 1.4 GHz flux from galaxies such as NGC 253, 
NGC 3628 and NGC 4631 by a factor ~ 2, where a relatively bright compact nuclear region 
of radio emission is surrounded by a fainter extended radio- emitting disk, i.e. cases where a 
single Gaussian flt the the flux distribution docs not work well. 

NGC 4945 is a southern hemisphere target and hence was not observed in the NVSS. 
For NGC 4945 we use the nuclear region 1.4 GHz flux and size quoted in Elmouttie et al. 
(1997). Their observations had a similar spatial resolution to that of the NVSS, making 
these values a reasonably good match to the values we use for the other galaxies. 

Note that the moderately-low spatial resolution of this radio data (FWHM = 45") will 
over-estimate the spatial size of the radio emission for the more compact starbursts (and 
hence under-estimate the SF intensity), especially for the more distant galaxies such as 
NGC 3079 and NGC 1482. 
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Fig. 1. — False color composite images of the central 5 kpc by 5 kpc of M82 (a), NGC 
1482 (b), NGC 253 (c), NGC 3628 (d), NGC 3079 (e), NGC 4945 (f), NGC 4631 (g), NGC 
891 (h) and NGC 6503 (i). Ha emission is shown in red, optical R-band in green and 
0.3-2.0 keV energy band soft X-ray emission is shown in blue. The X-ray emission is the 
point-source subtracted diffuse emission which has been adaptively smoothed to achieve a 
local S/N of 3, which tends to over-smooth structure in the X-ray emission. All images are 
shown on a square-root intensity scale. In contrast to Figs. 2 to 13 presented in Paper I, the 
same absolute intensity is not used, in order to show some of the fainter emission features. 
The galaxies are shown in order of (approximately) decreasing star formation rate per unit 
area. Note the close spatial similarities between the minor-axis-oriented Ha emission and 
the diffuse X-ray emission, in particular the limb-brightened outflow cones in NGC 1482, 
NGC 253, NGC 3079 and NGC 4945. 
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(a) 



•T-i 




Fig. 2. — As Fig. 1, except that panels (a) through (j) now show 20 kpc by 20 kpc boxes, 
while in panel (k) the region is 40 kpc by 40 kpc. The galaxies are M82 (a), NGC 1482 (b), 
NGC 253 (c), NGC 3628 (d), NGC 3079 (e and k), NGC 4945 (f), NGC 4631 (g), NGC 6503 
(h), NGC 891 (i), NGC 6503 (j), and NGC 3079 (k). 
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3. Summciry of results from Paper I 

In Paper I (Strickland et al. 2004) we presented a detailed study of the diffuse X- 
ray emission in a sample of fO approximately edge-on disk galaxies (7 starburst galaxies, 
3 normal spirals) that span the full range of star formation found in disk galaxies. With 
the arcsecond spatial resolution of the Chandra X-ray telescope we were able to separate 
the X-ray emission from point sources from the truly diffuse component. The high spatial 
resolution also allowed a meaningful comparison to be made of the spatial location of the 
X-ray emission with respect to multi-wavelength data (in particular to ground or space- 
based optical imaging). In addition to a detailed discussion of each galaxy, we presented a 
mini-atlas of soft and hard X-ray, Ha and R-band images of each of the 10 galaxies, shown 
at a common spatial and surface brightness scale to facilitate cross-comparison. Three color 
composite images of the diffuse soft X-ray emission (blue) with optical Ha emission (red) 
and R-band stellar continuum emission (green) are shown in Figs. 1 and 2. 

We presented a variety of quantitative measures of the spatial extent, spectral hardness, 
and shape of the diffuse X-ray emission, several of which can be applied to future samples 
of galaxies with lower S/N data. We found that the spectral properties of the diffuse X- 
ray emission show only weak, or some cases no, variation with increasing height (outside 
a heavily-absorbed, spectrally-harder region within z ~ 2 kpc of the disk mid-plane). The 
vertical decrease in surface brightness of the extra- planar emission (|2;| > 2 kpc) appears 
to better described by exponential (effective surface brightness scale heights are typically 
between 2-4 kpc), or Gaussian models, than the power law expected of a freely expanding 
fluid. 

For the eight galaxies with detections of extra- planar {i.e. halo- region) diffuse emission, 
we find that a common spectral model, comprising a two-temperature Mekal (Mewe et 
al. 1995) hot plasma model with an enhanced a-to-Fe element ratio, can simultaneously fit 
the ACIS X-ray halo-region spectra. The X-ray-derived metal abundances show super-Solar 
ratios of a-process elements (primarily oxygen) to iron. This is consistent with the origin of 
the X-ray emission being either (metal-enriched) merged SN ejecta, or from shocked ambient 
halo ISM (with moderate depletion of refractory elements onto dust). 

Although less luminous, the spatial and spectral properties of the thermal X-ray emission 
in the normal, non-starbursting, spiral galaxy NGC 891 are very similar to those of the 
starburst galaxies with superwinds. 

Our favored model for the origin of the extra-planar soft X-ray emission is that SN 
feedback in the disks of star-forming galaxies create, via blow out and venting of hot gas 
from the disk, tenuous exponential atmospheres of density scale height i7g ~ 4 - 8 kpc. 



-16- 



AGN-driven winds do not appear to be significant in this sample, in that there is no obvious 
correlation between the presence and luminosity of any AGN and the properties of the diffuse 
X-ray emission. The soft thermal X-ray emission observed in the halos of the starburst 
galaxies is either pre-existing halo medium (which has been swept-up and shock-heated 
by the starburst-driven superwinds) or from a small fraction (by volume) of the merged- 
SN-ejecta that has been compressed near the walls of the outflow (e.g. by a reverse shock 
propagating back into the outflowing wind). In either case the observed exponential X-ray 
surface brightness distributions are an inheritance from galactic fountain activity prior to the 
currently-observed starburst phase. This model is based on the qualitative 2-D morphology 
of the diffuse X-ray and optical Ha emission (in particular the fllamentary, occasionally 
limb-brightened morphology of both the X-ray and Ha emission) , as well as interpretation 
of the more- quantitative minor axis surface brightness and spectral hardness profiles. An 
implication of this model is that galactic-scale gaseous halos may be common around star- 
forming disk galaxies. Observing starburst galaxies, in which superwinds "light-up" these 
pre-existing halos, may present the best method for studying the gaseous halos of star- 
forming galaxies. 



Table 3. Summary and location of data values used in this paper 



Value Tabulated in Description Appears in 

(1) (2) (3) (4) 



Data values tabulated in this paper. 





Table 


1 


Total galactic IR luminosity, derived from IRAS fluxes. 


Figs 


3 and 4 


Lk 


Table 


1 


Total galactic K-band luminosity, derived from fluxes in the 2MASS Large Galaxy Atlas. 


Figs 


3 and 4 


Lb 


Table 


1 


Total galactic B-band luminosity. Not corrected for extinction. 


Figs 


3 and 4 


"rot 


Table 


1 


Inclination-corrected peak circular velocity. 


Figs 


5, 6 and 7 


/eo/Aoo 


Table 


1 


IRAS 60 to lOO^m flux ratio, an indicator of the intensity of star formation. 


Figs 


5 and 7 


Mtf 


Table 


1 


Galactic baryonic mass derived from the K-band Tully-Fishcr relation. 


Figs 


4 and 8 


-^^1.4GHz 


Table 


2 


Total galactic radio luminosity at 1.4 GHz, derived from the NVSS. 


Figs 


3 and 4 


/l.4GHz 


Table 


2 


Total galactic radio flux at 1.4 GHz, derived from the NVSS. 






^1.4GHz 


Table 


2 


Radio major-axis half light radius. 


Fig. 


6 


/l.4GHz/Si.4GHz 


Table 


2 


Radio-based estimate of the star formation rate intensity. 


Figs 


5 and 7 


/fir 


Table 


2 


Total galactic IRAS FIR flux. 






D25 


Table 


2 


Inclination corrected optical diameter at a surface brightness of 25 magnitudes per square arcsecond. 


Fig. 


6 


/FIR/-D25 


Tabic 


2 


Commonl\- used proxy of the mean galactic star formation rate per unit area. 


Figs 


5 and 7 


K-band t-q.s 


Table 


2 


K-band half light radius, from the 2MASS Large Galaxy Atlas. 


Fig. 


6 


M 


Table 


2 


Mean galactic mass surface density. 


Fig. 


6 




Table 


2 


Total galactic core-collapse supernovae rate, based on the total IR luminosity. 






-FSN,FIR,D25 


Table 


2 


One estimate of the mean supernovae rate per unit disk area, where i^sN,FIR,D25 — ^Sn/^25 


Fig. 


8 


fSN.FIR.Oi ,iGHz 


Table 


2 


Another estimate of the mean supernovae rate per unit disk area, where i^sN FIR T>2'^ ~ ^Sn/4''i 4GHz 


Fig. 


8 


Datii xidlLc,-, tiib.iliitc'il ill I'npri I, 


So. 5 


Table 


4 


Mean diffuse X-ray surface brightness between z = and the minor half light height zo.5 s-nd between —5 < r(fepc) < 5. 


Fig. 


7 


2^1e-9 


Table 


4 


Height along the minor axis at which the diffuse X-ray surface brightness reaches 10~^ photons s~^ cm~^ arcsec"^. 


Fig. 


8 


^0.95 


Table 


4 


Minor axis diffuse X-ray 95%-flux enclosed height. 


Fig. 


8 


'•0.5 


Table 


5 


Major-axis diffuse X-ray half light radius. 


Fig. 


6 


'■0.75 


Table 


5 


Major-axis diffuse X-ray 75%-flux enclosed radius. 


Fig. 


6 




Table 


6 


Minor-axis diffuse X-ray exponential scale height in the 0.3—1.0 keV energy band. 


Figs 


6 and 8 


< *T > 


Table 


8 


X-ray luminosity weighted mean temperature. 


Fig. 


6 


Ehalo 


Table 


8 


Mean diffuse X-ray surface brightness within the halo spectral extraction region. 


Fig. 


7 


/X.HALO 


Table 


9 


Halo region absorption-corrected X-ray flux (0.3—2.0 keV energy band). 


Fig. 


6 


i'X.HALO 


Table 


9 


Halo region absorption-corrected X-ray luminosity (0.3—2.0 keV energy band). 


Figs 


3 and 4 




Table 


9 


Disk region absorption-corrected X-ray luminosity (0.3—2.0 keV energy band). 






ix.NUCL 


Table 


9 


Nuclear region absorption-corrected X-ray luminosity (0.3—2.0 keV energy band). 












Sum of I^x,NUCL and ix,DISK- 


Figs 


3 and 4 



Note. — Column 2: Table in this paper, or paper I, in which the symbol in Column 1 is tabulated. Column 3: A short description of the mean of each symbol — see the appropriate 
Table for a more detailed explanation. Column 4: Figures in this paper in which the values are used. Some data values used in the figures, e.g. Z/x D+N' ^ot explicitly tabulated 
but can calculated from data that is tabulated. 
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4. Correlation analysis 

We shall now move onto a more-quantitative comparison of the properties of the diffuse 
X-ray with other properties of host galaxies. For convenience, definitions of the data values 
we use (from both Paper I and this paper), the table numbers from which they were extracted, 
and the figures in which they are used, are given in Table 3. In Paper I we considered the 
thermal emission from the nuclear region and the disk exterior to the nucleus separately. 
In this paper we will combine the X-ray luminosities of the nucleus (i^x,NUCL) and the disk 
(-^x,disk) to give a total X-ray luminosity Lx,d+n from hot gas with |^| < 2 kpc of the plane 
of the galaxy. For convenience we shall refer to this combination of Paper I's disk and nuclear 
regions as the disk in the rest of this paper. Halo region [i.e. \z\ > 2 kpc) luminosities are 
used as published in Paper I. 

Given the very significant absorption of the diffuse emission from NGC 4945 by the high 
Galactic foreground hydrogen, the low detection significance and atypical spectral hardness 
of the halo diffuse emission (see Paper I), we have excluded it from further analysis. In 
particular, it is excluded from Figs. 3 to Fig. 8. 

4.1. Star- formation or galaxy mass? The physical-drivers of extra-planeir 

X-ray emission 

The only previous study quantitatively comparing the diffuse X-ray emission^ from late- 
type, star-forming galaxies with their integrated luminosities at other wavelengths is the 
i?05^r PSPC-based study published in Read et al. (1997) and Read & Ponman (2001). 
They find that the estimated total diffuse emission X-ray luminosity per unit galaxy mass 
(estimated from Lx/ L-q) correlates well with Lb, an estimator of galaxy mass for normal 
spirals. For starburst galaxies Lx/ shows a better correlation with LpiR, i.e. with a proxy 
of star formation rate. 

There are a number of reasons while it is worthwhile to re-examine this issue, apart 
from the obvious advantage of using the higher quality data from an instrument with greater 
sensitivity, and vastly superior capabilities of separating point-source and diffuse emission. 
Determining the properties of extra-planar hot gas in normal galaxies has assumed a new 



^Significant effort has been expended in the past comparing the multi- wavelength luminosities of galaxies 
in an attempt to elucidate the physical causes of the integrated X-ray emission in these systems. We refer 
the interested reader to Fabbiano (1989), Shapley, Fabbiano & Eskridge (2001) and Fabbiano & Shapley 
(2002). 
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0.01 0.1 1 10 100 1000 10< 106 0.1 1 0.1 1 

L,, (10'" L,4Ghz(Lo) L^dO'^Lj L, (10'" L,,) 

Fig. 3. — Disk (including the nuclear region) and halo region diffuse emission 0.3-2.0 keV 
energy band X-ray luminosities (I/x,d+n and Lx.halO) open and filled circles respectively), 
plotted against the host galaxies' total IRAS IR luminosity (panel a), the 1.4 GHz luminosity 
(panel b), B-band optical luminosity (panel c) and NIR K-band luminosity (panel d). The 
dotted lines show the trend expected for a relationship of unit slope — they are not fits to 
the data. See § 4 for a definition of what is meant by disk and halo regions. 

significance, given its importance as a diagnostic of the accuracy of theories of galaxy for- 
mation and evolution. The normal galaxy sample of Read & Ponman (2001) is largely a 
collection of face-on systems, with the approximately Milky- Way-like galaxy NGC 891 (van 
der Kruit 1984; Wainscoat et al. 1987) being classed as a starburst on the basis of the 
reddening-correction-sensitive Lfir > 0.38 Lb criterion used by Read & Ponman as the 
definition of a starburst. It is only in high inclination systems where we can separate the 
disk and halo diffuse components, the latter being the most interesting component in testing 
both SN feedback models and cosmological accretion models. 

There are only three approximately edge-on normal spiral galaxies in our sample, the 
Milky Way-like galaxy NGC 891 and the two lower mass and lower luminosity spirals NGC 
6503 and NGC 4244. These three galaxies are too small a sample on their own to establish 
any hard conchisions about the properties of diffuse X-ray emission in normal spiral galaxies, 
but in combination with the starbursts wc can make a meaningful first assessment of how 
closely the properties of the diffuse X-ray emission in these systems resemble those of the 
edge-on starburst galaxies. Of course more X-ray observational study of edge-on spiral 
galaxies is required before we understand extra-planar hot gas in normal spiral galaxies as 
well as we understand the phenomenon in starbursts. 
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Fig. 4. — As Fig. 3, except all panels are luminosity per unit stellar mass. The disk and halo 
diffuse emission 0.3-2.0 keV energy band X-ray luminosities divided by total galactic stellar 
mass are shown and open and filled circles respectively. These values are plotted against IR 
luminosity per unit mass (panel a), radio luminosity per unit mass (panel b) and B-band 
optical luminosity per unit mass (panel c). The dotted lines show the trend expected for a 
relationship of unit slope. 



4.2. X-ray and multi-wavelength luminosities 

We find reasonable correlations between the soft X-ray luminosity of both the disk and 
halo diffuse X-ray emission, and the total galactic FIR and 1.4 GHz radio luminosities (see 
Fig. 3), which are well-known measures of the total star-formation rate. In particular, note 
that the slope of the correlations between Ljr or Liaguz and Lx is of order unity, as expected 
if the diffuse X-ray emission is caused by star-formation activity. 

Correlation between the X-ray luminosities and the optical B and NIR K-band lumi- 
nosities {i.e. proxies for total galactic stellar mass) is somewhat weaker, in particular in the 
case of the halo X-ray luminosity when only detections are considered. 

In Fig. 4 we plot all luminosities normahzed by the stellar mass of the host galaxy. 
Estimates of the stellar mass of each galaxy are based on the 2MASS K-band TuUy-Fisher 
relationship derived by Bell & de Jong 2001), and are listed in Table 1. As the K-band 
luminosity Lk and Mtf are not independent variables we have not included a plot of Lx/M^f 
against Lk/Mtf in this figure. Having taken out any dependence related to galaxy mass 
{i.e. more massive galaxies having more star-formation and more X-ray emission than small 
galaxies). Fig. 4 makes it clear that the diffuse X-ray luminosity per unit mass is directly 
proportional to star-formation rate per unit mass. 

In Fig. 5 we plot the ratio of the absorption- corrected X-ray flux from the halo to the 



- 21 - 



0.0001 



8xl0-s - 



6x10- 



o 

K 4xl( 



2x1 0-» h 




1 1 llllllj — 1 1 llll^ — 1 1 lllll^ 1 


1 1 lllll^ — 1 


1 iiiii^ 


: (a) 




- 


— 


X 


— 








T • 


X 


tj 






♦ : 


1 1 1 


1 1 


1 


0.0010.01 0.1 


1 


10 



0.0001 




fi.4 Ghz/^?.4 Ghz (mJy/arcsec2) 
0.0001 



100 150 200 250 

Vpoi (km/s) 



8x10-6 - 



6x10- 



o 

ffi 4X10-S 



I I I I I I I I I I I I I I I I I I I 
(c) 



4 T 



2x10-5 \- 5 



_L 



I ■ ■ ■ I 



200 



I 1 1 iiiiij 
(d)X 

I 150 P ^ 



0.2 0.4 0.6 0.8 1 

IRAS f,,/f,,. 



^ 100 - 

o 

< 

X - 
,x 50 - 





T 



uuiL 



uL 



iiL 



0.01 0.1 1 10 100 

fpiR/Dia (mJy/arcsec^) 



Fig. 5. — Panels a through d plot the ratio of the absorption-corrected 0.3-2.0 keV halo- 
region X-ray flux to the host galaxies IRAS FIR flux (a, b and c) or 1.4 GHz flux (d) against 
various proxies for the mean star formation rate per unit area in the host galaxy (panels 
a, c and d) and galaxy mass (panel b). Error bars show the la statistical errors due to 
uncertainties in the background-and-point source subtracted count rates alone - systematic 
uncertainties in /x,halo are probably a factor 2 or so. For NGC 6503 and NGC 4244 the 
3(7 upper limits are shown. For NGC 253 the extrapolated total halo-region flux is used. 
These results imply that the ratio of the halo X-ray luminosity to the star formation rate of 
the host galaxy is independent of the concentration or intensity of the star formation, or the 
mass of the host galaxy. 



host galaxy's FIR or radio flux, against estimators of the star-formation rate per unit disk 
area from radio and FIR observations, and the galaxy circular velocity. From panels a and b 
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of this figure, it is apparent that all galaxies with detected extra-planar X-ray emission have 
a similar fx,H alo/ Ifir ratio, of approximately 4 x 10~^ (±0.2 dex). The mechanical energy 
injection rate from supernovae and stellar winds in the starburst is Lsn ~ 0.01 e x LpiR, 
where e is the fraction of mechanical energy thermalized (Heckman et al. 1990). This 
implies that the halos of these starburst galaxies, along with the normal spiral galaxy NGC 
891, radiate an approximately fixed fraction ~ 0.004/e of the mechanical energy supplied by 
stellar feedback within the disk, in the 0.3-2.0 keV X-ray band. For any reahstic value of 
e, ie. 0.1 < e < 1, we re-derive the well-known result that superwinds radiate a negligible 
fraction of their energy budget in soft X-rays (we discuss this further in § 5.1.3). 

Perhaps the most interesting implication of Fig. 5 is that the X-ray luminosity per unit 
star formation rate is independent of the star formation intensity (star formation rate per 
unit disk area) , over many orders of magnitude variation in star formation rate per unit area. 
In other words, the efficiency of converting supernova mechanical power into X-ray emission 
appears to be independent of the original supernova rate per unit area or volume. 

It is also interesting that the normal spiral NGC 891 is so similar to the starbursts in 
this regard. Fig. 5 also makes it clear that the non-detections of diffuse emission in the halos 
of NGC 6503 and NGC 4244 are not particularly constraining - given their low FIR and 
radio brightness we would expect any diffuse X-ray emission generated by star-formation 
activity to be as faint or fainter than our existing 3a limits. 

4.3. The size of the X-ray emitting regions 

The major-axis X-ray half-hght radii ro.5 (tabulated in Table 5 of Paper I) provide a 
convenient measure of the radial size of the diffuse X-ray emission in and near the plane 
of each galaxy (as a large fraction of the diffuse X-ray emission comes from < 2 kpc of 
the plane). To reduce the influence of absorption on the measured sizes we plot ro.5 derived 
from the medium energy band (1.0-2.0 keV) against various simple measures of the radial 
extent of strong star-formation (^^i.4ghz)) the extent of the old stellar population (-D25 and the 
K-band ro.5), and the galactic mass (circular velocity frot) in Fig- 6a - d. The radial extent 
measure ro.5 correlates weakly with Oiaghz, D25 and K-band ro.5, and not at all with Vrot- 
The 75%-enclosed hght radius ro.75 shows a somewhat better correlation with D25 and the 
K-band ro.5. Quantitatively, the diffuse X-ray radial extent ro.5 is closest to the size of the 
non-thermal radio emission ^i.4ghz, as expected if the diffuse X-ray emission is generated by 
mechanical feedback from massive stars. The X-ray emission is more centrally concentrated 
than the radio emission, although this may due to the low spatial resolution of the NVSS 
data. 
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In contrast, the vertical extent of the diffuse X-ray emission in the halo appears in- 
dependent of ^i.4GHz- We have only shown the effective halo exponential scale height i^eff 
in Fig. 6e - h, but the other tracers of minor axis extent, such as ^0.95 or -Zie-g, produce 
very similar looking plots. Instead of correlating with local measures of where the energy is 
injected (such as ^i.4ghz), the distribution of X-ray-emitting gas at 2; > 2 kpc depends on 
the total size of the host galaxy, given the approximately linear relationship between Hcs 
and D25 and/or K-band ro.5. More quantitatively, Kendall's Tau statistic (Press et al. 1992) 
is only r = 0.21 for the correlation between ^i.4ghz and H^s {1.3a from zero correlation, 
equivalently probability of the correlation being spurious Ps — 0.19), whereas the stronger 
correlations mentioned above have r = 0.61 between D25 and H^s {3.Sa,Ps — 0.00017), 
r = 0.43 between the K-band ro.5 and ifeS {'^■7cr,ps — 0.0077) and r = 0.55 between v^ot and 
H^s (3.4(7, = 0.00076). 

Panel j of Fig. 6 plots Hes against the ratio of the luminosity-weighted mean halo 
temperature < kT > to the mean mass surface density in the disk = Mtf / D25 (note that 
the vertical component of the gravitational force is proportional to //). The correlation 
not particularly strong, in fact with a r = 0.46 (2.8cr for zero correlation, ps = 0.0045) it 
is weaker than the correlations between H^s and D25 or v^-ot- This suggests that the X-ray 
temperature < kT > and scale height H^s are not as directly related as would be the case 
for simple vertical hydrostatic equilibrium where H oc kT/ [i. 

We also investigated the relationship between iJcff and the ratio of v^^^/fi (panel k of 
Fig. 6). If the vertical support of the gas now visible in X-ray emission was originally due 
to semi-turbulent gas motions, then we might expect the vertical velocity dispersion to be 
proportional to the rotational velocity of the galaxy (as found for the stellar vertical velocity 
dispersion, e.g. van der Kruit & de Grijs 1999). Evidence in support of this scenario is 
lacking in the current data, as the correlation is weak, r = 0.31 {1.9a, Ps — 0.059). 
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Fig. 6. — A comparison of the characteristic physical dimensions of the diffuse X-ray emission 
to various proxies for the size of the star-burst region, the optical size of the host galaxy, 
and galaxy mass (the half-light radius of the 1.4 GHz radio emission Oi ^ghz, the inclination 
corrected D25 value from the RC3 catalog, the K-band half-light radius from the 2MASS 
Large Galaxy Atlas (Jarrct et al. 2003), and the circular rotation velocity Wrot)- In panels a 
to d these values are plotted against the 1.0-2.0 keV energy band major-axis diffuse x-ray 
half-light and 75%-enclosed light radii ro.5 (open symbols) and ro.75 (filled symbols), while 
in panels e through h they are plotted against the minor axis diffuse X-ray exponential 
scale height Hes in the 0.3-1.0 keV energy band. Data measured along the positive r or 
2;-direction are plotted as open triangles (ro.75 values are shown as filled circles), negative r 
or z as inverted triangles (ro.75 values are shown as filled squares). The dotted lines in panels 
a - c and e - g show the trend an exact one-to-one correspondence between the variables 
would produce. See discussion in § 4.3. In panels i through k we plot the vertical scale 
height against the logarithm of mean mass surface density (in Mq kpc~^), the ratio of the 
mean X-ray temperature to the mass density (units of keV kpc^ Mq~^), and the ratio of the 
square of the circular velocity to the mass density (km^ s~^ kpc^ Mq"^). 



-25- 



4.4. The brightness of the diffuse X-ray emission, compared to steir-formation 

intensities and galetxy masses 

The presence of correlations between the halo X-ray luminosity and total star-formation 
rate, and the vertical extent of the X-ray-emitting gas to the total size of the host galaxy, 
should automatically lead to correlations between the diffuse X-ray surface brightness and 
estimates of the star formation rate per unit area (star formation rate intensity, SFRI) in 
the host galaxy. 

This is indeed what we find, as shown in Fig. 7. In panels a through e we plot a variety 
of proxies for the SFRI, and in panel f the circular velocity, against the mean diffuse X- 
ray surface brightness within the half-light height E0.5. Panels g through 1 repeat this, this 
time using the halo region surface brightness Shalo obtained within the fixed physical-size 
apertures used in the halo-region spectral fitting. 

Although the half fight heights are typically only ~ 1 kpc or so, E0.5 is still a measure of 
the brightness of gas associated with the minor axis outfiow in the starburst galaxies of this 
sample. We have demonstrated in this and previous papers that the brightest diffuse X-ray 
emission is directly associated with optically-emitting gas for which unambiguous kinematic 
evidence of minor-axis outflow exists (e.g. , in Heckman et al. 1990; McKeith et al. 1995; 
Shopbell & Bland-Hawthorn 1998). In the normal spirals, NGC 891, NGC 6503 and NGC 
4244, there is no existing kinematic evidence for outflow, and hence one might expect S0.5 to 
be more of a measure of the brightness of the hot gas within the disk (as opposed to lying in 
projection within 2 kpc of the mid-plane), and also the mean brightness of unresolved point 
sources. 

Note that the correlation between S0.5, or Shalo, and the FIR SFRI fpiR/ D^^ (Fig- 7b 
and h, where the angular size Dfg is that of the old stellar population, and hence independent 
of the angular size of the main star forming region in the host galaxy) are better than the 
correlations between E0.5 or Ehalo and the NVSS 1.4 GHz SFRI fiAGujOl 4GHz (Fig- 7a and 
g, where ^i.4ghz is a direct, although low spatial resolution, measure of the size of regions in 
which SNe are occurring). This is due to the better correlation between the vertical extent 
of the extra-planar emission with the total size of the host galaxy, than with the radial size 
of the main star-formation regions (as discussed in § 4.3). 

The ordinate values in panels d, e, j and k of Fig. 7 are, to flrst order, star formation rate 
per unit galactic mass. These panels, along with panels f and 1 (which plot S0.5 and Ehalo 
against the circular velocity Vj-ot), show that the brightness of the diffuse X-ray emission 
does not correlate with galactic mass, and hence is genuinely caused by star formation. It 
is therefore not the case that the correlation in panels a through c, and g through h, is 
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spurious and arises simply because more massive galaxies, which can hold more hot gas in 
their potential wells, also have more star-formation in total. 

Again, it is intriguing that there is no obvious transition in X-ray surface brightness 
between the normal and starbursting spirals - the normal spirals follow the trends one could 
extrapolate on the basis of the starburst galaxies alone. The sample of non-starburst disk 
galaxies we have available to compare against the starbursts is very small, but no good 
reasons exist as yet to believe that these galaxies are peculiar. 

In summary, these results validate the basic assumption that the presence of extra- 
planar diffuse X-ray emission is driven by star-formation processes in the underlying disk, 
but that the large-scale structure of the host galaxy also plays a strong role in determining 
the properties and distribution of the resulting extra-planar emission. 

In Fig. 8 we show the observed vertical extent of the diffuse X-ray emission against 
estimates of the SFRI (recast in terms of supernova rate per unit area Fsn) and total galactic 
mass Mtf for the galaxies of our sample. We estimate the supernova rate based on the IRAS 
FIR fiux. Dividing this SN rate by the inclination corrected value gives -Fsn,fir,D25 i the 
mean SN rate per unit area of the entire disk. If instead we use the approximate diameter 
of the non-thermal radio emission 9ia we obtain FsN,FiR,ei 4ghz (See Tables 1 and 2 for the 
origin of these values). It is clear that there is no correlation between the vertical extent 
of the diffuse X-ray emission and the SN rate per unit area in the host galaxy, except that 
galaxies with lower SN rate per unit area than NGC 891 have no detected extra-planar X- 
ray emission. This is similar to other studies of extra-planar emission at radio and optical 
wavelengths. Dahlem, Lisenfeld & GoUa (1995) found no correlation between the vertical 
scale height of non-thermal radio emission and the star formation rate per unit area in a 
sample of edge-on disk galaxies (see also Dahlem et al. 2001). In a Ha imaging study of 
edge-on disk galaxies Miller & Veilleux (2003) also find no correlation between scale heights 
and the star formation rate per unit disk 

Fig. 6, in combination with Fig. 8, implies that the density scale height of the X-ray 
emitting material depends on the absolute size of the host galaxy, and not on the spatial 
distribution or intensity of star formation within it. Note that this is not a statement 
regarding the absolute physical extent of superwinds, which almost certainly is much greater 
than the region probed by current observations (and which will physically depend on the 
total amount of energy deposited by SNe and the density of the IGM). 

We were initially motivated to investigate the quantitative relationship between the 
extra-planar X-ray properties and the SFRI indicators LfirZ-DIs and /i.4ghz/^i,4ghz by opti- 
cal studies of the eDIG (see Dettmar 1992; Rand 1996; Hoopes et al. 1999; Rossa & Dettmar 
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2003, among others). These had demonstrated the presence of extra-planar optical emission 
correlated with Lfir/D^sj ^hc sense that galaxies with high values of Lpin/D^^ were more 
likely to show evidence of eDlG. Combined with the studies of Lehnert & Heckman (1995, 
1996); Dahlem et al. (2001) that demonstrate that IR-warm galaxies with /eo/Zioo > 0.4 
are extremely likely to have extra-planar emission, kinematic evidence of outflow and/or 
non-thermal radio halos, we expected a correlation between high star formation rate per 
unit area and the luminosity and extent of diffuse X-ray emission in disk galaxy halos. 

Nevertheless, in the X-ray band it appears that the correlation we have found between 
the surface brightness of the diffuse X-ray emission and star formation rate per unit area is 
the combination of two independent effects: 

1. A relatively uniform efficiency of converting mechanical energy to X-ray emission 
(§ 4.2), i.e. fx oc /fir, the exact physics of which is not fully understood. 

2. The vertical extent of the X-ray emission is directly related to the size of the host 
galaxy, Hes oc D25 or the K-band half-light radius (§ 4.3), possibly for the reasons 
discussed below. 

This does not mean that the SFRl is not an important parameter, in particular when 
considering the conditions necessary to drive gas into the halo of a galaxy (as we shall discuss 
below), but that other parameters play a significant role in shaping the emission we see. 

In terms of the interpretation of the extra-planar X-ray and Ha emission presented 
in Paper I (and summarized in § 3), item 2 above can be explained if disk galaxies have 
gas throughout their halos, with vertical density scale heights approximately similar to the 
optical diameter of the host galaxy. Depending on the exact model assumed the mass of gas 
in the halo of a Milky Way-like galaxy prior to any starburst is a few times 10'' to a few 
times 10^ Mq (see models 4 and 5 in Strickland et al. 2002). 

4.5. Limits placed on cosmological accretion models for the origin of hot gas 

in the halos of normal galcixies 

Without deep X-ray observations of a larger sample of edge-on normal spiral galaxies it 
is difficult to critically test cosmological accretion models of the type presented by Toft et 
al. (2002), and Sommer-Larsen et al. (2003). Excluding the starbursts from consideration 
for the time being being, we would expect halo 0.3-2.0 keV X-ray luminosities of order 
Lx = 4.60, 0.20 and 0.08 x lO^^ergs-^ for NGC 891, NGC 6503 and NGC 4244 respectively, 
based on the Toft et al. model. 
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These luminosities are based on converting the bolometric luminosity as a function of 
rotation velocity quoted in Toft et al. (2002) into Chandra energy band soft X-ray lumi- 
nosities. Here we have assumed the X-ray emission has a fixed temperature of 0.16 keV 
(the virial temperature for the NGC 891 halo kT^i-^ = nmnv^^^. The observed luminosity- 
weighted mean temperature obtained for NGC 891 is < kT >= 0.23 keV), a metal abun- 
dance Z = 0.3^0 (the conversion between the Lbol used by Toft et al. and 0.3-2.0 keV 
energy band Lx is largely independent of Z), and is described by an exponential surface 
brightness profile with a fixed surface brightness scale height of 2 kpc (as Toft et al. 2002 
do not provide any measure of how the characteristic size of the X-ray emission varies 
with galaxy rotation velocity). For such a distribution a fraction of the total lumi- 
nosity arises in gas at |2;| > 2 kpc. Toft et al. (2002) 's predicted bolometric luminosity 
(0.012-12.4 kcV energy band) of Lbol ~ 10^° (t;rot/230kms-i )5ergs-i (scatter ±0.5 dcx) 
thus converts to a total 0.3-2.0 keV X-ray luminosity {i.e. integrated over all heights z) of 
Lx ~ 1.4 X 10^^ (f rot/230 km s~^ )^ erg s~^ (for simplicity we ignore the fact that the halo 
temperature should vary with rotation velocity as approximately kT oc f^ot, as hence even 
lower X-ray fluxes would be expected in the 0.3-2.0 keV energy band for NGC 6503 and 
NGC 4244). 

The predicted accreted halo luminosities given in the first paragraph of this section 
should be compared the to the observed values of Lx,halo = 1-2 x 10^^ ergs~^ for N891 and 
3(7 upper limits of Lx.halo < 1-6 x lO^^ergs"^ and < 2.5 x 10^^ erg s^^ for NGC 6503 and 
NGC 4244 respectively (see Table 9 of Paper I) . This comparison indicates that the predicted 
accretion-model X-ray luminosity for a galaxy hke NGC 891 is close to that value observed, 
and that the NGC 4244 and NGC 6503 observations arc not deep enough to strongly constrain 
the accretion model. It is worth noting that it is far from clear what fraction of the gas 
accreted from the IGM by disk galaxies is heated to the halo virial temperature (see for 
example Katz et al. 2003 and references therein for a theoretical perspective). 

Although, based on luminosity alone, some fraction of the X-ray emission in the halo 
of NGC 891 might be due to gas accreted from the IGM, we find the reasoning presented 
in § 4.2, 4.3 and 4.4 more compelling as evidence that the extra-planar X-ray emission in 
NGC 891 is associated with, and dominated by, feedback from star-formation. This apparent 
commonality of normal spirals with starbursts, based on the properties of their extra-planar 
emission, might be an artifact of the small number statistics. Deep X-ray observations 
of other edge-on normal spirals, spanning a broad range of star-formation rates, are clearly 
necessary before firmer conclusions can be drawn. Nevertheless, a robust finding of a smooth 
transition in the halo gas properties from normal to starburst galaxies would have major 
implications for understanding the large-scale structure of the ISM and stellar feedback 
processes in galaxies. 
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Fig. 7. — Diffuse X-ray surface briglitness plotted against various proxies for tlie intensity 
of star formation (star formation rate per unit disk area from /i.4ghz/6'i 4ghz' /firZ-DIs or 
log/eo/Zioo)) total star-formation rate per unit galaxy mass (Lir/Lb or somewhat more 
accurately Lir/Mtf), and galaxy mass (as Mtf oc t'roD- Panels a - f plot the effective 
X-ray surface brightness E0.5 within the half-light height zq.s, which is largely a measure of 
the brightness near or within the disk. Open symbols correspond to the 0.3-1.0 keV energy 
band, filled symbols to the 1.0-2.0 keV energy band. Triangles and squares represent data 
measured along the positive- 2; direction, inverted triangles and circles data measured along 
the negative- ;2 direction. Panels g - 1 plot the mean halo region surface brightness in the 
0.3-1.0 keV band (open circles) and the 1.0-2.0 keV band (filled circles). The six-sided stars 
are cases where values from positive and negative z overlap. See discussion in § 4.4. 
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Fig. 8. — The vertical extent of the diffuse X-ray emission as a function of the mean SN rate 
per unit area Fg^ or circular velocity v^ot of each galaxy. See § 5.1.4 for a description of how 
-fsN.FiR,D25 Ei'iid -^SN,FiR,6»i 4ghz Were calculated. Panels a - c plot the fitted surface brightness 
scale height for z > 2 kpc, data from positive z shown as open triangles and negative z as 
inverted triangles. The horizontal dashed line corresponds to a gas density scale height of 
2 kpc. Panels d - g plot the isophotal minor axis size Zie^g (positive z data points shown 
as open triangles, negative z as inverted triangles) and 95% enclosed minor axis light height 
-^0.95 (positive z data points shown as filled circles, negative z as filled squares). Where no 
measurement can be made a negative value is shown. 
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5. Compcirison the theoretical models of mechanical energy feedback from 

massive stars 

5.1. Disk blow-out energy requirements 

In a seminal paper, Mac Low & McCray (1988) considered, both analytically and 
numerically, the amount of mechanical energy that supernovae and stellar winds need to 
supply for a superbubble to blow-out of the stratified atmosphere of a galactic disk. Blow out 
occurs if the dense shell of swept-up gas marking the outer boundary of a superbubble ceases 
to decelerate and begins to accelerate down the density gradient. At this point Rayleigh- 
Taylor instabilities set in and rapidly fragment the superbubble shell, leading to the venting 
of the hot high pressure gas from the interior of the bubble out into the surrounding medium. 
As this venting of hot gas (the gas blowing out) occurs when the outer shell fragments and 
breaks up the term "breakout" is often used synonymously with blow out. This theory, 
with minor alterations, has been widely used since that time, in particular with respect 
to assessing the the prospects for metal ejection from galaxies (e.g. Koo & McKee 1992; 
Martin 1996; Mac Low 1996; Tenorio-Tagle 1996; Silich & Tenorio-Tagle 1998; Mac Low & 
Ferrara 1999; Ferrara & Tolstoy 2000; Ferrara, Pettini & Shchekinov 2000; Silich & Tenorio- 
Tagle 2001). Note that blow out of a superbubble from the disk ISM of a galaxy does not 
automatically imply that gas escapes into the IGM — more mechanical work must be done 
by the vented hot gasses to move any gaseous halo medium. Simulations of superwinds show 
that these hot gasses blow new bubbles in the halo, sweeping the ambient halo gas up into 
new shells (e.g. Tomisaka & Bregman 1993; Suchkov et al. 1994; Strickland & Stevens 
2000). We shaU address this in § 5.2. 

Mac Low & McCray demonstrate that whether or not a bubble blows/breaks out 
of a given atmosphere depends only on the mechanical energy input rate, i.e. mechanical 
luminosity, driving the bubble, and the pressure and scale height of the ambient ISM the 
bubble is expanding into — for convenience we re-derive this result below. For such an 
important theoretical construct, there has been surprisingly little in the way of observational 
testing and calibration. Heckman et al. (1990) demonstrated that the energy injection rates 
in starburst galaxies exceeded this critical blow-out luminosity (hereafter Lcrit)j but is is 
not clear from that work (and other comparisons to starbursting galaxies e.g. Martin 1996) 
what the lowest mechanical energy injection rate is that leads to blow-out in real galaxies. 

In this section we shall use the observationally determined presence of hot gas in the 
halos of our sample galaxies as empirical evidence for blow-out from the disk, and compare 
the estimated rates of SN energy return in disks of these galaxies to theoretical values of 
LcRiT- At present the observational results do not strongly constrain the theory. When larger 
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samples of edge-on spirals have had deep X-ray observations made of them the promise of 
this method should be realized. 



5.1.1. The critical mechanical luminosity for disk blow-out 

The theoretical critical mechanical luminosity required for blow-out is based on the 
following concepts. If, in a stratified disk with a vertical gas density distribution — 
Po exp{—z/H), an expanding superbubble or SN remnant has an expansion velocity Vbub 
(the shell expansion velocity) greater than the sound speed Cg = {'^P/pY^'^ in the ambient 
ISM by the time it reaches a vertical height of order the scale height H of the disk, then the 
dense shell will begin to accelerate. At this time the bubble or remnant shell will fragment 
due to Rayleigh- Taylor instabilities. The hot, high pressure, plasma interior to the shell will 
vent out through the fragmented shell into the surrounding medium, i.e. the halo for blow 
out from the disk of a galaxy. The dense fragments of the superbubble shell will be carried 
out along with the hot gas, in fact the ram pressure of the hot gas streaming past them will 
accelerate the shell fragments to higher velocities than the original shell expansion velocity 
^^bub (see e.g. Chevaher & Clegg 1985; Heckman et al. 1990, 2000). 

If instead I'bub ^ Cs, the bubble's expansion will continue to decelerate and eventually 
coast to a stop without blowing out. Radiative cooling will eventually lead to depressuriza- 
tion of the bubble interior, followed by the collapse of the bubble, in addition to damage 
caused by galactic differential rotation and random ISM motions. 

A more formal analysis identifies I'bub > Cg at z = 3H as the critical requirement for 
blow-out in a single- phase exponential atmosphere (see Ferrara & Tolstoy 2000), numerically 
verified by simulations showing blowout occurring once the bubble reaches a height of ~ 3H 
(Mac Low & McCray 1988; Mac Low, McCray & Norman 1989) 

For a bubble driven by a constant mechanical luminosity Lw into a medium of uniform 
and constant ambient density po (Castor, McCray & Weaver 1975; Weaver et al. 1977), the 
radius rbub and shell expansion velocity Vbub as a function of time t are 

r^.^^aL'^'p-'/'t'/', (1) 

and 

^;bub = 0.6aL^^'po'/'r2/^ (2) 

where a = (125/1547r)^/^. Rearranging these equations gives us the expansion velocity at 
the point that the bubble radius is 3H: 

^3H = 0.6 X 3-^/^ ^„ a'/' Po'^'H-'/^ (3) 



-33- 



where ^„ is a numerically-determined correction factor to account for the vertical decrease 
in density leading to bubble expansion velocities greater than in the constant density case. 
Using a code that numerically integrates the appropriate equations of motion (Strickland & 
Stevens 1999), we find = 2.5l 

The resulting minimum mechanical luminosity required to blow out of this atmosphere, 
-^CRiT, is 

LcRiT = (1547r/3) 7^/^ ^"^ P,"' pT H\ (4) 

where Pq is the thermal pressure in the ambient ISM. Note that this luminosity is a factor 
~ 3 lower than the critical luminosity derived by Mac Low & McCray (1988), i.e. blow-out 
from the disk is slightly easier. Converting to units more appropriate for the disk of a normal 
spiral, 

T / D \ 3/2 

-t^CRIT 



r--3n„-(^) (^), (5) 



]^Q38 gj^g g 

where tiq = pQ/fimu and = lO^^'^gcm"^. More energy is needed to blow out of a 

starburst region, where measured ISM pressures reach P/k ~ 10^ K cm^^ (McCarthy et al. 
1987; Heckman et al. 1990; Lehnert & Heckman 1996). 



5.1.2. Disk blow out due to individual stars clusters, or through cooperation between 

multiple clusters? 

With respect to blow out from the disk of normal star forming galaxies is it commonly 
assumed that each superbubble is powered by a single star cluster. In contrast, it has 
long been recognized that superwinds from galaxies with nuclear starbursts are powered 
by the collective energy return from the massive stars within the entire starburst region, 
typically a few hundred parsecs in size. This is a well established observational result that 
can be demonstrated in various ways, e.g. the superwind pressure profiles in Heckman et 
al. (1990) and Lehnert & Heckman (1996) demonstrate that mass and energy are being 
injected relatively uniformly over a region equal in size to the starburst region. Another 
example is that the edges of best-resolved superwind outflows, in galaxies such as M82 and 
NGC 253, matching up well to the edges of the nuclear starburst regions (Strickland & 
Stevens 2000; Strickland et al. 2000). This is by no means surprising, as star clusters 



^Ferrara & Tolstoy (2000) present a purely analytical derivation of fan including the effect of the density 
decline, but we find this method over-predicts vsu by a factor 3.7 over our numerical result (even when errors 
in their equations 24 and 25 arc corrected for). We speculate that this discrepancy is due to their neglect of 
the inertia of the cold shell in their calculation of the shock velocity. 
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tend themselves to cluster spatially just as massive stars tend to form within clusters, with 
starbursts representing the extreme of the phenomenon (one can describe the M82 starburst 
as ~ 100 R136's within a region only 2-3 times larger than 30 Doradus, see O'Connell et 
al. 1995). Furthermore, Meurer et al. (1995) demonstrate that only a fraction ~ 20% of 
the massive stars within nearby starburst regions are associated with very massive clusters, 
there being a substantial population of massive stars in starburst regions not associated with 
obvious clusters (see also Tremonti et al. 2001; Chandar et al. 2003). 

The question arises whether superbubblcs can be cooperatively driven by massive stars 
in multiple stellar clusters in situations less extreme than those in nuclear starbursts, i.e. in 
disk-wide starbursts or even in normal star forming galaxies? This process need not be as 
efficient or as common as in the nuclear starbursts, but it might occur often enough to be 
significant in terms of the disk/halo interstellar interaction. There are a number of reasons 
why it is worthwhile to consider such a scenario. Firstly, individual clusters that are massive 
enough to drive a superbubble capable of blowing out of the disk of a normal spiral galaxy 
on their own are relatively rare, as we show in § 5.1.2. If only one or two such clusters exist 
in a galaxy at any given time it would be difficult to sustain a galaxy-wide fountain-driven 
hot halo (assuming that these are not extremely peculiar spiral galaxies). Secondly multi- 
dimensional simulations of the ISM with energy return by SNe and stellar clusters show 
that bubbles do merge, and once merged can potentially become powerful enough to escape 
the disk (de Avillez 2000). A third issue, related to the first, is that the normal spiral 
galaxy NGC 891 has a very similar ratio of halo- region X-ray luminosity to FIR luminosity 
{i.e. total star formation rate) to the starbursts with superwinds, implying a similar total 
efficiency in converting the initial supernova mechanical energy into soft X-ray emission. If 
cooperative action does not occur in normal spirals and only the most massive cluster drive 
bubbles that blow out of the disk, then a smaller fraction of the massive stars forming in 
NGC 891 would be responsible for the hot gas its the halo. To achieve a /x.halo//fir ratio 
similar to that in the starbursts would then require fine tuning the physical processes in the 
halo of NGC 891 to be more efficient at radiating soft X-rays, despite the other similarities 
to the extra-planar X-ray emission from the starbursts. 

A final reason to consider blow out due to collective superbubble driving is that it 
naturally leads to an analytical formulation in which there is a critical star formation rate (or 
supernova rate) per unit area necessary for blow out. From the stand point of observational 
studies of galaxies beyond the local group it is much easier to estimate a mean star formation 
per unit area than to measure the properties of individual stellar clusters. Indeed, individual 
star clusters are likely be heavily obscured in the edge-on galaxies where the presence of 
extra-planar hot gas is empirical evidence for blow out from the disk. Furthermore estimates 
of the star formation rates in these galaxies may only be available using low spatial resolution 
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FIR observations. As we shall discuss later, this formulation remains useful when considering 
blow out from the halo even if cooperative action does not apply to superbubble blow out 
from the disk. 

Energy return from non-interacting star clusters 

Let us now consider what the mass is of an OB association that contains enough massive 
stars enough to power a superbubble capable of blowing out of a disk of a normal spiral 
galaxy, assuming that only a single cluster can power a specific superbubble. The critical 
energy given in Eqn. 5 can be recast in terms of the number of massive stars formed in 
the cluster, and from that in terms of star cluster mass. Assuming all stars in a cluster 
form instantaneously, then the period over which energy injection occurs corresponds to the 
lifetime of the least massive star to undergo core collapse, i.e. stars with initial masses Mj ~ 
8Mq, which is ~ 30 - 40 Myr. The energy injection rate over this period is approximately 
constant, see e.g. Leitherer & Heckman (1995). The time needed for a bubble with Lw = 
-^CRiT to expand to the point it breaks out of the disk is similar to this time period, validating 
the assumption of continous energy injection in the formulae used to derive Lcrit- We 
assume that each SNe returns e x 10^^ erg of mechanical energy to the ISM, where the 
thermalization efficiency e is the fraction of the initial kinetic energy per supernova that can 
be used to drive a bubble, i.e. the energy fraction not radiated away prior to each young 
remnant overlapping with other SN remnants. We further assume a Salpeter IMF with 
upper and lower mass limits of 100 and 1 Mq respectively, which gives one massive star per 
~ 53 Mq of stars formed. Where we to assume a more realistic IMF extending down to the 
hydrogen burning limit (e.g. that of Kroupa, Tout & Gilmorc 1993) we would get less SNe 
per mass of stars formed, i. e. the cluster masses required for disk blow out would be larger 
and hence blow out rarer. 

The Milky Way has a gaseous scale height of ~ 0.4 - 0.5 kpc in the neutral ISM (Dickey 
& Lockman 1990; Sembach & Danks 1994), and between 0.5 - 1.5 kpc for the free electrons 
tracing the warm ionized component (Reynolds 1989; Nordgren, Cordes & Terzian 1992). 
The mean number density and thermal pressure in normal star-forming spiral like the Milky 
Way or NGC 891 are Uq ~ Icm-^ and P/k ~ 10^ Kcm-^ (Ferriere 2001; Wolfire et 
al. 2003). Assuming H = 0.5 kpc, Uq = Icm"^ and P/k = 10^ Kcm"^ then Lcrit = 
8.6 X 10^^ erg s~^ in a typical region of a normal spiral galaxy. This correspond to the energy 
return from 109e~^ SNe, or a cluster mass of Mcl ~ 5.8 x 10^e~^ Mq. A canonical value of the 
supernova thermalization efficiency under normal interstellar conditions is e = 0.1 (Thornton 
et al. 1998, note that e is expected to be higher in regions of higher star formation rate per 
unit volume such as star clusters e.g. Larson 1974; Canto et al. 2000). These results are 
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consistent with those of Ferrara et al. (2000), who also considered blow out from the disk 
assuming that superbubbles can only be blown by a single star cluster. For galaxies with 
star formation rates equivalent to the normal spiral galaxies we are discussing they derived 
a critical number of SNe for blow out similar to, or larger than, the number we derive above. 

Star clusters as massive as Mql are not common in normal, non-starbursting, galaxies. 
By way of comparison consider the masses of the most spectacular and massive star clusters 
in the Local Group. The Arches cluster in the Galactic center, although containing ~ 5% 
of all the Wolf-Rayet stars in the Galaxy, has a mass < 7 x IC^Mq (Figer et al. 2002). 
The R136 cluster in the 30 Doradus nebula of the LMC has an estimated initial mass of 
> 1.7 X 10^ M0 (Malamuth & Heap 1994), or ~ 3 x 10^ M©, upper limit 1.5 x 10^ 
(Brandl et al. 1996). The entire Orion Nebula Cluster, has a total gravitating mass of 
4.5 X 10^ Mq (Hillenbrand & Hartmann 1998). From this it is clear that the most massive, 
most populous, star clusters are capable of driving superbubbles powerful enough to blow 
out of normal spiral galaxy disks on their own, unless e is of order unity. 

How common are such clusters, and what fraction of massive stars form in these massive 
clusters? The number of massive stars per star cluster appears to follow a power law such 
that the number of of clusters with N massive stars is proportional to (see Oey et al. 
2003, and references therein). Wc will treat this as a continuous distribution with limits 
NuiN = 1 and A^max for computational convenience, as the results from the exact solution 
(see Oey ct al. 2003) are essentially identical in the regime we consider. The fraction of star 
clusters containing N or more massive stars is thus />jv = {N~^ ~^max)/(-'- ~^max)- 
generously assume that even normal galaxies can potentially form star clusters with masses 
up to 10^ Mq, then A^max ~ 1.9 x 10^. Assuming a lower A^max will reduce />Ar. For cluster 
thermalization efficiencies of 0.1, 0.3 and 1.0 the required number of SN for disk blow out are 
1090, 327 and 109 respectively, with corresponding cluster masses of Mcl = 6 x 10^, 1.8 x 10^ 
and 6 x 10^ Mq. For these three efficiencies the fraction of star clusters powerful enough to 
blow out of the disk of a normal spiral galaxy are />7v = 8.7 x 10~^, 3.0 x 10~^ and 9.1 x 10~^ 
respectively. In other words it is rare to find a star cluster that is capable of driving, on its 
own, a superbubble that can blow out of the disk of a normal spiral galaxy. 

Collective energy return over galactic scales 

Motivated by the preceding arguments we explore the requirements for blow out from the 
disk of normal spirals in a scenario at the opposite extreme from isolated, non-interacting 
star cluster, model commonly adopted. We shall assume that separate star clusters can 
cooperatively blow a single bubble, just as starburst regions drive a single superbubble and 
then superwind. This gives rise to a critical star formation rate (or supernova rate) per 
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unit area necessary for blow out, which is a number that is easier to empirically constrain in 
edge-on disk galaxies (e.g. see Table 2) than the masses of specific clusters. 

Mac Low & McCray (1988) find that at blowout the radius of the super bubble within 
the disk is ~ H. We therefore assume that some significant fraction T of all SNe due to 
massive star formation within an area ttH^ of the disk can contribute energetically to creating 
a blowout, i.e. . 

Ly^^nH^Fs^eTEs^, (6) 

where Fsn is the SN rate per unit disk area, and Eq-n — 10^^ erg is the kinetic energy return 
per SN event. 

Note that hydrodynamical simulations demonstrate that the expansion of a wind-blown 
bubble or superbubble is relatively insensitive to the spatial distribution of the energy in- 
jection within it (e.g. Strickland & Stevens 1998, 2000), for the reasons discussed by Mac 
Low & McCray (1988). We can thus continue to use the standard expansion laws despite 
having a more randomly distributed source of mechanical power. 

Physically e is the thermalization efficiency within individual clusters where radiative 
losses are small (e.g. Canto et al. 2000), while T represents the effective efficiency of com- 
bining the energy output from the multiple clusters within the nH^ area of disk that can 
possibly contribute to a bubble just about to blow out. Theoretically, both e and T are func- 
tions of the star formation rate, as supernovae alter the porosity of the ISM, and hence alter 
radiative energy losses and the effectiveness of collective action (Larson 1974). Relatively 
little theoretical work has been done with the aim of estimating e for reahstic cluster condi- 
tions, and this theoretical work has not been tested against observation in any meaningful 
way. Although simulations of the ISM that include interacting superbubbles do exist (e.g. 
de Avillez 2000), these tend not to discuss the energetic efficiency with which separate bub- 
bles merge. One situation somewhat similar to that we envisage can be found in a numerical 
simulation of a primordial galaxy by Mori, Ferrara & Madau (2002). They simulated bubbles 
blown by 90 star clusters spread over a region ~ 1 kpc in radius and measured the fraction 
of the total mechanical energy injected by SNe escaping into the IGM to be ~ 30 - 50%, 
depending on the initial conditions used. Their simulation effectively assumed e = 1, and 
their galaxy is not exactly comparable to a present day spiral disk (e.g. it lacks a halo), but 
these results imply T > 0.3 - 0.5 (in our model eT is the fraction of available SN mechanical 
power released within the disk that is available to drive the final bubble that blows out, but 
less than this fraction of energy will reach the halo due to radiative energy losses that would 
occur even in the absence of superbubble merging and collective action). One would ex- 
pected T to be of this order in normal disk galaxies. Physically the primary radiative energy 
losses in two merging bubbles would come from the colliding shells, as these are dense and 
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hence will be efficient radiators if shocked in the coUision. However, only ~ 20% of the total 
mechanical energy originally fed into a superbubble is thought to be contained in the shell, 
the majority of the energy (~ 45%) residing in the hot interior and ~ 35% being lost during 
shell formation (Mac Low & McCray 1988). The hot gas in merging bubbles is unlikely to 
suffer significant additional radiative losses, given its low density. For starbursts, where the 
evidence for effective collective action is best, T may be of order unity, as the higher SN 
rate per unit volume should lead to higher ISM porosity and hence efficient thermalization 
(e.g. Larson 1974). In other words the net efficiency eT throughout the entire starburst 
region might be comparable to the thermalization efficiency immediately within a dense star 
cluster (e.g. Canto et al. 2000). By way of comparison the SN rate per unit volume in the 
center of M82 is approximately two orders of magnitude higher than the SN rate per unit 
volume in the simulations of Mori et al. (2002) that we discuss above. 

Combining Eqn. 6 with Eqn. 4, we obtain the critical SN per unit disk area for blow 
out from the disk, Fsn,crit, to be 

FsN,CRiT = (154/3) r''^ p,^l^ Pl'^ (e T E^^yK (7) 

Note that this is no longer dependent on the scale height. Large gaseous scale heights do 
not inhibit blow out as the bubble can draw mechanical energy from a correspondingly larger 
region of disk. Scaling to more convenient units, 

FsN,CRiT(Myr-^kpc-^) = 3.45 x n''^' (^^^ ^ , (8) 



or alternatively 



VlO^A; J \eT J \10^ J 



/ ^ \ 1/2 



FsN,CRiT(Myr-kpc-) = 3.45 x ( ( ^) ( ^) . (9) 



For the purposes of comparing observation with theory Fsn should ideally be calculated 
on a location-by-location basis within a given galaxy, based on the local supernova rate 
T^-sN, or star formation rate, summed within an area ir H^, where H is the appropriate gas 
density scale height for the region we are interested in assessing blow out from (be it the thin 
molecular disk, the more extended Reynolds layer, or the halo itself) — i.e. some knowledge 
of the scale height is still required. 



5.1.3. Constraints on feedback efficiencies from the observed luminosities 



Both the standard bubble blow out theory presented in § 5.1.1 and the slightly modified 
version above depend strongly on the efficiency of mechanical energy feedback through terms 
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such as e and T. Theoretical work suggests these efficiencies depend on the local star 
formation rate through the porosity of the ISM (Larson 1974), but the exact values are hard 
to quantify and there is little agreement between different groups (e.g. contrast Strickland 
(2002) with Recchi, Matteucci & D'Ercole (2002)). 

Adopting a different approach, we can crudely constrain the net mechanical energy 
feedback efficiency by relating the rate of mechanical energy injection into the halo to the 
observed X-ray luminosities. From the relationship between the FIR luminosity and SN 
rate given in Heckman et al. (1990), i?sN = 0.2Lfir/10^^ Lq, we obtain that the rate of 
mechanical energy injection into the halo is 

i^w,HALO~ 0.0167 eTLpiR, (10) 

where we ignore energy losses due to superbubble shell formation and any superbubble 
merging that occurs (see § 5.1.2). Including these radiative losses might reduce Lw.halo by 
a factor ~ 2, and lead to a similar increase in the estimated efficiencies, if analytical theory 
and numerical models of winds are right (see e.g. Strickland & Stevens 2000). 

Direct observation of the volume-filling high-pressure gas that drive superwinds would 
be the best method of measuring the supernova thermalization efficiency. In the absence 
of this we are forced to estimate e from a theoretically-determined relationship between 
the X-ray luminosity and the mechanical power actually driving the wind. The fraction 
^XRAD = -^x,HALo/-^w.HALO dcpcuds ou the physical process responsible for the soft X-ray 
emission, e.g. non-radiative wind shocks, or possibly conductive interfaces between hot and 
cold gas. For X-ray emission from conductive interfaces (Weaver et al. 1977; Chu & Mac 
Low 1990), which is the mechanism that has been most extensively explored to date, Axrad 
depends only weakly on the ambient density and the value of the thermal conductivity (Lx 
scales almost exactly in proportion to Lw, and increases weakly for higher ambient density 
or thermal conductivity). We calculated Axrad for conditions appropriate for superwinds 
using the code described in Strickland & Stevens (1999), which numerically solves the 
Weaver et al. (1977) model and uses the Mewe et al. (1995) hot plasma code to calculate 
the X-ray emission. For Lw,halo in the range 10^° to lO'^^ergs"^ , gas densities in the range 
10~^ to 1 cm~^ , ages in the range 5 to 30 Myr and values of the thermal conductivity in 
units of the Spitzer conductivity 0.1 to 10, we find a mean ^xrad = 0.013 (with a scatter 
±0.6 dex). In our hydrodynamical simulations of superwinds (Strickland & Stevens 2000), 
which do not include thermal conduction, we find a similar fraction of the mechanical energy 
emerging in the soft X-ray band (again with a significant scatter, depending on the exact 
initial conditions assumed). 

For the time being we shall assume that Axrad ~ 0.01 is a reasonable estimate of the 
fraction of the wind mechanical power that is radiated in the soft X-ray band. Combining 
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this with Equation 10 we obtain 

^x,HALO - 1.67 X 10-^ (Axrad/0.01) eT Lpm- (11) 

The typical /x,halo / /fir ratio we observe is ~ 4 x 10"^, with the range being 2.2 - 8.0 x 10~^, 
which this imphes eT ~ 0.24 with a range of 0.14 - 0.48. The normal spiral galaxy NGC 
891 has one of the lowest /x,halo//fir ratios in the sample, with /x,halo//fir = 2.2 x 10~^, 
and thus a lower than average value for eT ~ 0.14 (NGC 891 is the lowest if we exclude the 
somewhat doubtful detection of halo emission from NGC 4945, although it is only a factor 
~ 2 lower than the mean /x,halo//fir ratio), a lower than average value for eT ~ 0.14. 

5.1.4- Comparison to observations 

If, as seems most probable, the origin of the diffuse extra-planar X-ray emission around 
disk galaxies is due to SN feedback, then empirically assessing whether significant blow out 
from the disk has occurred is simple. It effectively reduces to identifying which galaxies 
have diffuse X-ray emission at z >2 kpc {i.e. several disk scale heights). With a sample of 
edge-on galaxies, some of which have extra-planar hot gas and other which do not, and with 
reasonable estimates of the star formation or SN rate per unit area for each galaxy we can 
then empirically determine (a) whether there is a critical SN rate per unit area (Fsn,crit) 
above which blow out from the disk occurs, and if so (b) the value of this critical SN rate 
per unit area. 

In principle we could compare the empirically derived value of Fsn,crit with the theo- 
retical value (e.g. Equ. 8), and hence test whether the theory is correct or not. In practice it 
is currently impossible to use this method to meaningfully test superbubble blow out theory 
due to a mixture of theoretical and observational uncertainties. 

The primary theoretical uncertainties are in the mechanical energy thermalization effi- 
ciency for individual star clusters e and efficiency of cooperation between multiple clusters T. 
The largest uncertainties in these parameters exist the for star forming conditions in. normal 
disk galaxies, i.e. galaxies that probably span the critical threshold in star formation rate 
per unit area for disk blow blow out. 

Observational difficulties include uncertainties in estimating the star formation or su- 
pernova rates, and more significantly the uncertainties in the effective area and hence the 
mean SN rate per unit area. 

Currently it is only practical to estimate the total supernova rate T^-sn of each of 
these galaxies (or other edge-on disk galaxies that will be observed with Chandra or XMM- 
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Newton) . This prevents us from calculating Fsn in the location- by-location manner described 
in § 5.1.2^. Nor do we have exact measures of the scale height of the disk ISM, and Fsn 
depends sensitively on the assumed scale height. Considering the issue of disk blow out, the 
ISM scale height is likely to be of order a few hundred parsecs to a kiloparsec or two. As 
this is similar to the magnitude of ^^i.4ghz (see Table 2), one could argue that it is reasonable 
to use 7^sn/(46'^ 4qjj2) as a plausible Fsn to use when assessing disk blow out (this value is 
denoted as Fsn,fir,6»i 4ghz Table 2 and in Fig. 8) . In contrast the D25 diameter (typically 
~ 10 - 20 kpc) is commonly used when estimated mean galactic star formation rates per unit 
area in optically-based studies of extra-planar warm ionized gas (e.g. Rand 1996; Hoopes 
et al. 1999; CoUins et al. 2000; Rossa & Dettmar 2003; Miller & Veilleux 2003). This leads 
to estimated mean SN rates per unit area orders of magnitude lower (this value is denoted 
as FsN,FiR,D26 in Table 2 and in Fig. 8). 

From Fig. 8 it appears that all galaxies with Fsn.fir.^i 4ghz ~ 2000 SN Myr~^ kpc~^ 
have diffuse extra-planar X-ray emission, indicative of disk blow out. Let us assume that the 
divide between disk blow out and superbubble confinement occurs at or near this value of 
Fsn- As the galaxy with the lowest Fsn value that also has detected extra-planar emission is 
the normal spiral galaxy NGC 891, let us assume that no ~ lcm~^ , P/k ~ lO^Kcm"^ and 
eT ~ 0.14. Thus equation 8 would then predict disk blow out for Fsn ^ 25 SN Myr~^ kpc~^, 
a factor ~ 80 lower than the observationally-estimated value, apparently a strong discrepancy 
between blow out theory and observations. However, if we use Fsn,fir,D25 instead, then 
galaxies with Fsn,fir,D26 ^ 40 SN Myr~^ kpc"^ appear to have disk blow out, apparently in 
good agreement with blow out theory. 

In summary, although it is relatively easy to empirically determine that a galaxy has 
blown hot gas out of the disk into the halo, a quantitative comparison to theory is much 
more difficult. With the present uncertainties it is not clear whether the energy requirements 
for superbubble blow out that the standard theory predicts (with or without allowing for 
collective energy return from multiple clusters) are accurate, even to within an order of 
magnitude. More theoretical and observational study is required before we can meaningfully 
test superbubble blow out theory: 



^We note that Ha imaging could be used to evaluate i^sN on a location-by-location basis within the same 
galaxy, although this method has typically been used on to calculate mean galactic SFRI values and blow-out 
thresholds (e.g. Rossa & Dettmar 2003). Ha emission in the disk is typically much more radially extended 
than the regions responsible for driving the outflow (the exception being M82, where star-formation is now 
confined solely to the nuclear region), i.e. star formation in some regions of the disk is occurring at rates per 
unit area below the critical value. This is graphically illustrated in Fig. 1 and 2. Of course, Ha measurements 
suffer the weakness that they systematically underestimate the true SFRI in the regions with the highest 
star formation rates due to the associated higher levels of extinction. 
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1. X-ray observations must be deep enough to rule out significant amount of hot gas in 
the halos of low SF rate galaxies, so that we can be certain that blow out from the 
disk has not occurred (see § 4.2). 

2. Observational estimates of star formation or supernova rates per unit area should ide- 
ally use spatially resolved estimates of the star formation rate, and more importantly, 
observed disk gaseous scale heights (see § 5.1.2). 

3. We need more theoretical study of the energetic efficiencies of mechanical feedback 
in clusters and in the realistic multi-phase ISM, in order to understand how e and T 
depend on the star formation intensity and local ISM properties (see § 5.1.3). 

5.2. Ejection of metal-enriched gas into the IGM by superwinds 

Much of the preceding discussion as revolved around the blow out of superbubbles from 
the disk of normal star forming galaxies. We now wish to briefly turn to the issue of ejection 
of hot, metal-enriched, gas into the IGM from starburst-driven superwinds. 

5.2.1. The plausibility of halo blow-out 

Material blowing out from the disk does not automatically escape into the IGM, as it 
must work against the obstacle provided by gas in the halo before being able to reach the 
IGM, as mentioned in § 1. 

Based on the observed X-ray surface brightness profiles of the extra-planar emission of 
the galaxies within this sample, it is appropriate to represent the original halo gas distribution 
as an exponential with density scale height i7g ~ 4 - 8 kpc, i.e. twice the observed surface 
brightness scale height. Note that optical studies of the extra-planar emission from warm 
ionized gas in normal and some starburst galaxies also indicate exponential gas distributions 
with very similar scale heights (Hoopes et al. 1999; CoUins et al. 2000; Miller & Veilleux 
2003). 

For such a stratified halo gas distribution the standard superbubble blow out energy 
requirements can be applied (this time with respect to halo blow out) to any bubble or wind 
that manages to blow out of the disk. Simulations of superwinds clearly show that the halo 
medium is swept-up into a new shell surrounding the hot gas vented from the disk (Suchkov 
et al. 1994; Strickland &; Stevens 2000). We will assume that superbubble blow out theory 
is correct for the remainder of this section. Note that this directly implies that the critical 
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region of space around a galaxy that can prevent escape into the IGM is within ~ 3Hg (~ 12 
- 24 kpc) of the mid-plane, a region that can be probed by existing observational methods. 
We need not concern ourselves about the density and pressure of gas at larger, ~ 100 kpc, 
distances from the host galaxy^. 

We can now apply the arguments of § 5.1 to halo blow out. The proportionality between 
the halo region X-ray flux and the IRAS FIR flux (and hence the total galactic SF rate) 
strongly suggests that the halo can draw on mechanical energy from a significant fraction 
of the entire star-forming disk. If the SN rate per unit area in the disk is sufficient to blow 
out from the disk, then the collective blow-out model can be applied to the hot gas now in 
the halo. The appropriate scale height to use is that of the halo gas density distribution 
Hg, either directly measured from the X-ray data or estimated using empirical results of 
§ 4.3. There we showed that Hg oc D25 (or the K-band half hght radius ro.5). Thus the 
appropriate mean supernova rate per unit area Fsn to calculate for assessing halo blow out 
is proportional to TZsn/D^5 = -^sn,fir,D26 ■ 

The only significant differences from the case of assessing disk blow out are then as 
follows: 

1. The appropriate density and pressure are that of the halo, evaluated at the galactic 
mid-plane (e.g. see ShuU & Slavin 1994). For typical galactic environments now and 
at high redshift, the density and pressure of any bound, non-transient, halo medium 
will be less than those in disk. 

2. The additional energy losses during disk blow must be accounted for. In the star- 
bursts, a good case can be made that radiative losses are negligible, based both on 
the observational lack of significant X-ray, O VI and optical emission (Heckman et al. 
1990, 2001; Hoopes et al. 2003). In the high resolution hydrodynamical simulation of 
Strickland & Stevens (2000) > 50% of the mechanical energy injected within the disk 
is transported into the halo. 

3. When considering blow out from the halo it is appropriate to use a SN rate aver- 
aged over the characteristic size of the host galaxy, not just the main star forming 
region. For comparison to Eqn. 7 a good estimate of Fsn would be 7^sN/('?r-ffg) ~ 
'^sn/ (47rrQ 5 j^_^,3^^jj) or 7^sN/(0.047rD25), given the correlations shown in Fig. 6. 



^This is true unless the galaxy happens to be in a high density environment such as a cluster or compact 
galaxy group, where the pressure and density of the inter-cluster or inter-group medium may be significant, 
in which case the arguments of Silich & Tenorio-Tagle (1998); Murakami & Babul (1999); Silich & Tenorio- 
Tagle (2001) apply. 
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Making full use of this to assess blow out from the halos of galaxies requires that we know 
the critical SN rate per unit area for disk blow out. Nevertheless, with an observationally- 
calibrated theory of this sort it will be possible to assess whether gas can blow out into 
the IGM. What this analysis does indicate is that extended gaseous halos around normal 
star-forming galaxies (excluding galaxies in dense gaseous environments such as clusters) do 
not automatically prevent outflows from reaching the IGM. The total power requirements 
for blowing out of a gaseous halo will more somewhat more exacting than for simply blowing 
gas out of the disk, but this is not surprising, nor a significant problem given the SN rates in 
nuclear starbursts. Nevertheless the most important point we wish readers to extract from 
this discussion is that the critical region of space around a galaxy that can prevent escape 
into the IGM is most probably within ~ 3Hg (~ 12 - 24 kpc) of the mid-plane, which is a 
region well probed by existing observational methods. 



5.2.2. Requirements for escape into the IGM 

Unambiguous observational proof of ejection of material into the IGM by outfiows re- 
quires kinematic evidence of outward motion ai v > v^sc at heights z > 3Hg above the 
mid-plane of the host galaxy (where v^sc is the escape velocity at that location, and not 
the escape velocity from the center of the halo, and Hg is the density scale height of the 
gaseous halo)'^ ^. This is observationally challenging, especially given that the most interest- 
ing, i.e. metal-enriched, gas phases are very tenuous, and very hot and highly ionized. Such 
observations, even of the denser phases, are unlikely to be accomplished for many years. 
Nevertheless, what the theoretically-and-observationally motivated argument in the preced- 
ing section demonstrates is that it is likely that starburst galaxies, even ones as massive as 



''Note that v, and preferably the summed thermal and kinetic energy, must be assessed on a phase-by- 
phasc basis. The cool shell velocity (as used by, for example Ferrara & Tolstoy 2000: Fcrrara ct al. 2000), 
or the temperature of the X-ray emitting plasma (e.g. Wang et al. 1995; Martin 1999), are not particularly 
meaningful tools for assessing whether gas will be retained. The energy per particle in the hot phases is much 
greater than that of the cool shell material, and it is likely that the kinetic energy of the X-ray-emitting hot 
gas can exceed its thermal energy several times (Strickland & Stevens 2000). 

^Wc wish to emphasize that gas motions at velocities less than the local escape velocity do not necessarily 
mean that the gas will be retained by the host galaxy. The gaseous phases in superwinds for which kinematics 
are available are entrained into the wind, either by Kelvin- Helmholtz instabilities along the wall of the outflow, 
or as fragments of the original superbubble shell that fragmented at disk blow out. Their motions are not 
ballistic, as they are embedded in the higher velocity fluid of merged supernova cjccta. Following blow out 
from the disk their kinematics are no longer determined by the standard Weaver et al. (1977) wind-blown 
bubble model. 



-45 - 



M ~ 10^^ Mq, and possibly some normal spiral galaxies, do eject material into the IGM. This 
is consistent with the effective yields and metal-loss scenario presented by Garnett (2002) 
and Tremonti ct al. (in preparation). Put another way, we regard escape as an eminently 
plausible fate for at least some fraction of the material in superwinds, as currently there are 
no convincing theoretically or ohservationally-motivated ways to absolutely prevent escape. 

6. Summeiry 

Making use of the detailed measurements of the diffuse X-ray emitting gas in and around 
the sample of seven starburst and three normal edge-on spiral galaxies presented in Paper 
I, we have investigated how the properties of this hot gas correlate with the size, mass, star 
formation rate and star formation intensity in the host galaxies. Although not a large sample, 
these galaxies do span the full range of star formation activity found in spiral galaxies. The 
three normal spirals are too small a sample on their own to establish any hard conclusions 
about the properties of diffuse X-ray emission in normal spiral galaxies, but in combination 
with the starbursts we can make a meaningful initial assessment of how closely the properties 
of the diffuse X-ray emission in these systems resemble those of the edge-on starburst galaxies. 

We demonstrate that the extra- planar {i.e. halo-region) diffuse X-ray emission in the 
starburst galaxies is ultimately driven by star formation activity within the disk, i.e. through 
starburst-driven superwinds. Accretion of gas from the IGM, or AGN-driven winds (see 
Paper I), do not appear to be significant in this sample. 

Considering both the starburst and normal spiral galaxies, larger galaxies show more 
radially and vertically extended diffuse X-ray emission, but beyond this correlation we find 
no evidence that galactic mass plays any part in determining the diffuse hot phase in the disks 
and halos of these galaxies. The extent of the diffuse X-ray emission parallel to the plane 
of the host galaxy correlates well with estimates of the extent of star formation within the 
disk, further evidence that the soft thermal X-ray emission is a result of mechanical energy 
feedback from massive stars. The vertical extent of the diffuse X-ray emission correlates best 
with the absolute size of the host galaxy, as if the emitting gas fills the halo. Exponential 
scale heights are similar to those found in the extra-planar warm ionized gas from optical 
for samples of star forming disk galaxies. 

The luminosity of the diffuse X-ray emission, both within the disk and halo, is de- 
termined primarily by the rate of mechanical energy injection due to SNe in each galaxy 
(which is correlated to the star- formation rate) . The surface brightness of the extra-planar 
diffuse X-ray emission correlates with the star formation rate per unit area in the underlying 
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disk, but this correlation is a combination of the previously mentioned luminosity and size 
correlations. 

The properties of the diffuse X-ray emission in the normal spiral galaxies is very similar 
to that found in the starburst superwind galaxies. In fact, one could successfully predict the 
diffuse X-ray properties of all three normal galaxies in this sample based on extrapolating 
the X-ray properties of the starburst galaxies to the lower star formation intensities seen in 
the normal galaxies. The close similarity in the halo-region diffuse X-ray emission in the 
starburst (with superwinds) and the Milky- Way-like normal spiral NGC 891 (most probably 
having a star-formation-fed galactic fountain) is intriguing, and deserves more study. Is it 
simply due to similar, micro-physical, emission mechanisms, or does the similarity extend 
to the macroscopic hydrodynamics at work over multi-kpc scales? The current sample of 
edge-on normal spiral galaxies observed by Chandra or XMM-Newton is small, so it is still 
possible that it is purely a coincidence that the three normal spirals we have studied lie on 
trends extrapolated from the starbursts. With the larger sample of observations available 
with next few years we should finally be able to typify the properties of diffuse thermal X-ray 
emission in normal spirals in the way we have with starburst galaxies. 

We also consider what the observed extra-planar diffuse X-ray emission can tell us 
about the both the efficiency of massive star mechanical feedback on galactic scales, and the 
accuracy of standard theoretical models of superbubble blow out from disk galaxies. 

We find the ratio of the halo-region diffuse X-ray flux to the host galaxy's total FIR flux 
is typically /x,halo//fir ~ 4 x 10~^, to within a factor two, for both the starburst galaxies 
and NGC 891. This ratio is independent of the star formation rate per unit area or the mass 
of the host galaxy. 

It has long been clear that starburst-driven superwinds are powered by the collective 
mechanical energy input from massive stars within the entire starburst region. In contrast 
superbubble blow out from the disk ISM of normal star forming galaxies is usually considered 
in the context of bubbles being powered by a single star cluster, in which only very rare, 
very massive, star clusters arc powerful enough to drive a bubble that can blow out of a disk 
(Ferrara & Tolstoy 2000; Ferrara et al. 2000). We argue that there are a variety of good 
reasons to consider some level of cooperative or collective action between multiple clusters in 
normal star forming galaxies, i.e. when considering blow out of gas from the disk a significant 
fraction of all SN energy with a disk area ~ nH^ can be used to blow a superbubble. We 
show that a minor generahzation of the Mac Low & McCray (1988) superbubble model to 
allow for such collective action directly predicts a critical star formation rate per unit area 
for superbubble blow out from the disk, a formulation particularly useful for comparison to 
observations of edge-on star forming galaxies. 
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The standard theory of superbubble blow out developed by Mac Low & McCray is 
widely used in astronomy, largely in relation to mass and metal loss from galaxies, but has 
had little in the way of observational testing. The presence of hot gas in the halo of a 
star forming galaxy is direct evidence for blow out of superbubbles from the disk of the 
host galaxy (in cases where the emission can be shown to be related to SF activity in the 
underlying disk). Further observational study of edge-on normal spiral galaxies in the ways 
presented in this paper, and in Paper I, will allow empirical constraints to be made on the 
critical star formation rate per unit disk area necessary to blow hot gas out of the disk into 
the halo. In principle these empirical measurements could be used to test superbubble blow 
out theory. We discuss and emphasize several observational and theoretical issues that must 
be overcome before superbubble blow out theory can be practically tested. 

Finally, we apply this collective-action blow out theory to superwinds blowing out from 
the extended gaseous halos of their host galaxies, as our and other workers observational 
results strongly suggest the gaseous halos of star forming galaxies are well represented as 
exponential atmospheres. It is important to note that blow out from the halo into the 
IGM is energetically more difficult that blow out from the disk into the halo. The existing 
uncertainties in the accuracy of blow out theory prevent robust quantitative assessment of 
whether the metal enriched gas powering superwinds escapes into the IGM. Nevertheless, 
we argue that is is still probable that the hot, metal-enriched gas in superwinds does escape 
into the IGM, even for galaxies as massive as M ~ 10^° - 10^^ M© in normal environments 
{i.e. not in dense galaxies groups or clusters). Irrespective of the quantitative uncertainties 
associated with disk and halo blow out, we wish to emphasize that the distribution of gas in 
the halos of star forming galaxies observed in the optical and in the soft X-ray band implies 
that the crucial spatial region around a galaxy that controls whether gas will escape into 
the IGM is not the outer halo ~ 100 kpc from the host galaxy. Instead it is the inner few 
halo scale heights, within ~ 20 kpc of the galaxy plane, a region well probed by existing 
observational methods. 

Within the next few years it will be possible to calibrate theories of disk and halo blow 
out and mechanical energy feedback against observations of local star-forming galaxies, and 
thus make meaningful inferences about the enrichment of the IGM by outflows, and the 
influence of massive star feedback on cosmological scales. 

Over the several years this project has been in the making, we have have been fortunate 
to beneflt from the insightful comments and questions of many astronomers, too numerous to 
mention individually, to whom we extend our thanks. We would also like to thank S. Hameed 
for providing us with R-band and Ha images of NGC 1482. We thank the anoymous referee 
for comments that lead to a clearer presentation of our arguments. 
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